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ABSTRACT 


(Distribution  Limitation  Statement  A) 

A  computer  program  has  been  developed  for  use  in  determining 
the  dynamic  response  of  an  aircraft  to  runway  roughness  during  takeoff 
and  constant  speed  taxi.  The  mathematical  model  has  been  programmed 
in  Fortran  for  a  CDC  6600  digital  computer.  A  typical  takeoff 
simulation  requires  less  than  200  seconds  of  computer  time  and  less 
than  77,000  octal  storage  locations.  The  output  from  the  program  is 
in  two  formats,  a  digital  listing  and  a  Calcomp-plotted  time  history. 
The  plotted  output  is  very  useful  in  evaluating  results. 
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SECTION  I 


INTRODUCTION 

This  volume  contains  a  detailed  presentation  of  the  fora  and  use  of 
the  program  "TAXI",  which  will  compute  the  dynamic  response  of  an 
aircraft  to  a  runway  profile.  This  program  consists  of  one  main 
program  and  four  subroutines.  Input  to  the  program  consists  of  airplane 
data  on  punched  cards  and  a  runway  profile  on  magnetic  tape.  Its  output 
includes  a  listing  of  ten  aircraft  parameters  at  given  time  intervals 
and  a  time  history  Calcomp  plot  of  center  of  gravity  and  pilot  station 
vertical  accelerations  and  the  runway  profile  traversed  by  the  nose  gear 
of  the  aircraft. 

This  program  will  enable  the  user  to  simulate  any  aircraft 
traversing  a  . unway  profile  with  a  minimum  of  modification  to  the 
computer  code.  The  phrase  "conventional  landing  gear  aircraft"  refers 
to  aircraft  with  one  set  of  main  landing  gear  which  are  nonarticulated 
and  of  the  single-acting  type  such  as  the  C-141,  KC-135  and  B-52. 

Several  nonconventional  landing  gear  aircraft  include  the  F-4,  F-lll 
and  C-5A.  All  conventional  landing  gear  aircraft  can  be  simulated  with 
only  a  change  in  the  aircraft  input  data.  In  contrast,  the  F-4,  vhich 
is  not  a  conventional  landing  gear  aircraft,  requires  a  change  in 
aircraft  data  and  the  substitution  of  two  subroutines  in  the  basic 
computer  program  with  those  designated  for  use  for  the  F-4.  Simulation 
of  the  F-lll  also  requires  a  change  in  aircraft  data  and  substitution 
of  two  subroutines  in  the  basic  computer  program  with  those  designated 
for  the  F-lll.  Due  to  the  complexities  of  modifying  the  basic  program 
by  substitution  of  subroutines  for  the  C-5A  aircraft,  which  has  two 


sets  of  main  gear  that  are  double-acting,  a  completely  separate  deck 
is  used. 

In  order  to  facilitate  the  use  of  the  program  TAXI,  this  volume 
also  contains  a  detailed  description  of  the  program  in  its  basic  form 
and  the  modifications  of  the  basic  source  deck  required  for  the  F-4, 
F-lll  and  C-5A  simulations.  First,  the  form  of  the  aircraft  input  data 
cards  and  runway  profile  tape  are  described  along  with  the  general 
setup  of  the  program  source  deck.  Next,  program  flow  charts  are  shown. 
Finally,  the  output  listing  and  plot  are  discussed.  A  complete  listing 
of  the  program  and  Fortran  symbols  are  contained  in  appendices  1  and  2. 


SECTION  II 


PROGRAM  OPERATION  PROCEDURES 

General  Procedures 

Before  the  program  TAXI  can  h«  used  for  a  runway  roughness  analysis, 
the  following  procedures  must  be  followed. 

First,  adequate  computer  core  memory  space  must  be  specified.  For  all 
simulations  a  core  memory  of  77000  octal  words  is  required. 

Second,  sufficient  computer  time  must  also  be  specified.  Central  pro¬ 
cessor  (CP)  time  is  a  function  of  the  integration  step  size,  number  of  flexible 
modes  included  in  the  analysis,  and  the  type  of  simulation,  either  a  taxi  or 
takeoff.  In  order  to  estimate  this  CP  time,  the  following  examples  can  be 
used.  Using  the  recommended  integration  step  size  of  0.001  sec.,  a  takeoff 
with  10  modes  of  vibration  included  requires  CP  time  of  2.5  sec.  for  every 
second  of  actual  simulation  time.  For  a  taxi  with  10  modes  included,  CP  time 
is  2.0  sec.  for  every  second  of  actual  simulation  time.  Including  fewer  vib¬ 
ration  modes  and  a  larger  step  size,  the  result  will  be  less  CP  time  required 
while  more  modes  and  a  smaller  integration  step  size  will  require  more  time. 
These  estimates  are  for  a  CDC  6600  computer.  Other  machines  may  require  more 
or  less  time. 

Third,  a  set  of  data  cards  for  the  aircraft  being  simulated  must  be 
punched  -»nd  a  magnetic  cape  on  which  the  runway  profile  is  stored  must  be 
obtained.  The  formats  for  the  aircraft  data  cards  and  magnetic  tape  are  con¬ 
tained  in  Section  II  of  this  Volume. 
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Fourth,  the  modifications  of  the  basic  source  deck  must  be  made  if 
an  F-4,  or  F-lll  simulation  is  to  be  run.  If  a  C-5A  aircraft  is  to  be 
simulated,  the  source  deck  designated  for  this  particular  aircraft  must 
be  used.  Instructions  for  the  modification  of  the  basic  source  deck 
for  an  F-4  or  F-lll  simulation  are  contained  in  Section  II  of  this 
volume. 

Finally,  the  entire  deck  including  control  cards,  source  deck 
and  aircraft  data  cards  must  be  assembled  as  shown  in  this  section. 


are  necessary  for  an  F-4  simulation.  All  other  cards  remain  the  same. 

Table  III  contains  the  changes  to  various  cards  in  Table  I  for  an 
F-lll  simulation.  All  other  data  cards  remain  the  same. 

Similarly,  Table  IV  contains  the  changes  to  certain  data  cards  in 
Table  I  required  for  a  C-SA  simulation.  Again,  all  other  data  cards 
remain  unchanged. 


The  aircraft  input  data  for  the  C-141,  KC-135,  B-52,  F-lil,  F-4 
and  C-5A  aircraft  are  contained  in  Volume  I  of  this  report.  A  sample 
aircraft  input  data  deck  for  the  KC-135  aircraft  is  shown  in  Table  VI. 


5 


TABLE  I 


AIRCRAFT  DATA  FOR  CONVENTIONAL  AIRCRAFT 


Section  1  (cards  1-5) 

-  General  Airplane  Data 

Card 

Column 

Format 

Variable 

Name 

Definition 

Card  1 

1-80 

8A10 

PLANE 

Airplane  Being  Simulated  and  Gross  Weight 

Card  2 

1-10 

F10.1 

W 

Vehicle  Weight  (lbs) 

11-20 

F10.1 

A 

Distance  Main  Gear  to  CG  (in) 

21-30 

F10.1 

B 

Distance  Nose  Gear  to  CG  (in) 

31-42 

F12.0 

MMI 

Mass  Moment  of  Inertia  (lb  in  sec^) 

Card  3 

1-10 

F10.2 

PSARM 

Distance  of  Pilot  Scation  to  CG  (in) 

11-20 

F10.2 

TAILRM 

Distance  of  Tail  Station  to  CG  (in) 

Card  4 

1-10 

F10.2 

SPEED 

Initial  Velocity  of  Airplane  (ft/sec) 

11-20 

F10.2 

THRUST 

Total  Airplane  Thrust  (lbs) 

21-30 

F10.2 

TAKOFF 

Airplane  Rotation  Speed  (ft/sec) 

Card  5 

1-10 

F10.4 

CL 

Lift  Coefficient 
'~‘>*g  Area  (ft2) 

11-20 

FIG.  4 

AREA 

21-30 

F10.4 

CD 

Drag  Coefficient 

Section 

2  (cards  6-11) 

-  Main  and  Nose  Gear 

Card  6 


1-10 

F10.2 

WM 

Unsprung  Weight  cf  Each  Main  Gear  (lbs) 

11-20 

F10.2 

WN 

Unsprung  Weight  of  Nose  Gear  (lbs) 

21-30 

F10.2 

SXM 

Number  of  Main  Gear  Struts 

31-40 

F10.2 

SXN 

Number  of  Nose  Gear  Struts 
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TABLE  1  (Con't) 


Card 

Co lumn 

Format 

Variable 

Name 

Definition 

Card  7 

1-10 

FI  0.5 

AHN 

Hydraulic  Piston  Area  Nose  (in2) 

11-20 

F10.5 

AAN 

Pneumatic  Piston  Area  Nose  (in^) 

21-30 

F10.5 

AHM 

Hydraulic  Piston  Area  Main  (in2) 

31-40 

F10.5 

AAM 

Pneumatic  Piston  Area  Main  (in2) 

Card  8 

1-10 

F10. 5 

PA0N 

Nose  Strut  Preload  Pressure  (Ibs/in^) 

11-20 

F10.  5 

PAOM 

Main  Strut  Preload  Pressure  (lbs/in2) 

21-30 

F10.5 

VON 

Fully  Extended  Nose  Strut  Air  Volume  (in^) 

31-40 

F10.5 

VOM 

Fully  Extended  Main  Strut  Air  Volume  (in^) 

41-50 

F10.5 

0  AM 

Orifice  Area  Main  (in2) 

51-60 

F10. 5 

OAN 

Orifice  Area  Nose  (in2) 

Card  9 

1-10 

F10.3 

SLH 

Distance  front  Axle  to  CG  Waterline 

Main  Gear  Strut  Unloaded  (in) 

11-20 

F10.3 

SLN 

Distance  from  Axle  to  CG  Waterline 

Nose  Gear  Strut  Unloaded  (in) 

Card  10 

1-10 

F10.1 

TSMI 

Main  Tire  Spring  Constant  Per  Strut  (Ibs/in) 

11-20 

F10.1 

TSNI 

Nose  Tire  Spring  Constant  Per  Strut  (ibs/in) 

Card  11 

1-10 

F10.5 

DX 

Integration  Step  Size 

Card  12 

1-5 

15 

IFF  LOT 

0  Plot 

1  No  Plot 


Section  3  (cards  13-16)-Metering  Pin  Description 


Card  13 

1-5 

15 

NSCN 

Number  of  Slope  Changes  Nose  < 

7 


TABLE  I  (Con't) 


Card 

Variable 

Co]  uran 

Format 

Ncane 

Definition 

*Card  14A, 

14  B, . . . . 

1-10 

F10.3 

STROKh 

Stroke  Corresponding  to  Metering  Pin  Diameter, 
Nose  Gear 

11-20 

F10.3 

PINDN 

Metering  Pin  Diameter,  Nose  Gear  (in) 

Card  15 

1-5 

15 

NSCM 

Number  of  Slope  Changes  Main  Gear 

*Card  16A, 

1 6B .  i 

1-10 

F10.3 

STROKM 

Stroke  Corresponding  to  Metering  Pin  Diameter, 
Nose  Gear 

11-20 

F10.3 

PINDM 

Metering  Pin  Diameter,  Main  Gear  (in) 

Section  4  (cards  17-19)-Flexibllity  Data 


Card  17 

1-5  II 

**Card  18A,  18B, 


Number  of  Flexible  Modes 


1-10 

F10.3 

SIMA1N 

Mode  Shape  Deflection  Main  Gear 

11-20 

F10.3 

SINOSE 

Mode  Shape  Deflection  Nose  Gear 

21-30 

F10.3 

SICG 

Mode  Shape  Deflection  CG 

31-40 

F10.3 

SITAIL 

Mode  Shape  Deflection  Tail  Station 

41-50 

F10.3 

SIPS 

Mode  Shape  Deflection  Pilot  Station 

**Card  19A, 

19B, .... 

1-15 

F15.2 

GM 

Generalized  Mass  (lbs  sec^/in) 

16-25 

FI  0.3 

OMEGA 

Modal  Frequency  (rad/sec) 

*Ore  card  is  required  for  each  stroke-metering  pin  combination  read  into  the 
program. 

**0ne  card  is  requited  for  each  flexible  mode. 
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TABLE  II 


INPUT  DATA  CHANGES  FOR  THE  F-4 
Card  Variable 

Column  Format  Name  Definition 


Card  13 


1-5 

15 

NSCN 

Number  of  area  changes  in  nose  gear  metering 
tube 

*Card 

14A,  l^B,.... 

1-10 

F10.3 

STROKN 

Stroke  corresponding  to  orifice  area,  nose 

11-20 

F10.3 

PINDN 

gear 

Net  orifice  area  at  STROKN,  nose  gear  (in^) 

*Card 

11 

1-5 

15 

NS  CM 

Number  of  area  changes  in  main  gear  metering 
tube 

*Card 

16A,  16B,,... 

1-10 

F10.3 

STROKM 

Stroke  corresponding  to  orifice  area,  main 

11-20 

F10.3 

PINDM 

gear 

Net  orifice  area  at  STROKM,  main  gear  (inz) 

*One  card  is  required  for  each  stroke-metering  pin  combination  read  into  the 
program 
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TABLE  III 


INPUT  DATA  CHANGES  FOR  THE  F-lll 


Card  Variable 

Column  Format  Name  Definition 


Card  13 


1-5  15  NSCN 

*Card  14A,  14 B . 


Number  of  area  changes  in  nose  gear  fluted 
metering  pin 


1-10  F10.3  STROKN 

11-20  F10.3  PINDN 


Card  lb 


1-5 


15 


NS  CM 


*Card  36A,  16Bt.... 


Stroke  corresponding  to  orifice  area,  nose  gear 
Net  orifice  area  at  STROKN,  Nose  Gear  (in2) 


Number  of  area  changes  in  main  gear  fluted 
metering  pin 


1-10 

11-20 


F10.3  STROKM  Stroke  corresponsing  to  orifice  area,  main  gear 

F10.3  PINDM  Net  orifice  area  at  STROKM,  main  gear  (in2) 


*  One  card  is  required  for  each  stroke  -  metering  pin  combination  read  into  the 
program. 
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TABLE  IV 


INPUT  DATA  CHANGES  FOR  THE  C-5 


> 


Card 

Column 

Format 

Variable 

Name 

Definition 

Card  2 

1-10 

Fl.0.1 

W 

Vehicle  Weight  (lbs) 

11-20 

F10.1 

A 

Distance  from  Rear  Main  to  CG  (in) 

21-30 

F10.1 

B 

Distance  from  Nose  Gear  to  CG  (in) 

31-40 

F10.1 

r* 

Distance  from  Front  Main  to  CG  (in) 

41-52 

F12.0 

MM  I 

Mass  Moment  of  Inertia  (lb-in-sec2) 

Card  13 

1-5 

*Card  14A, 

15 

14B  a •  ♦  •  • 

NSCN 

Number  of  Area  Changes  in  Nose  Gear 
Metering  Tube 

1-10 

F10.3 

STROKN 

Stroke  Corresponding  to  Orifice  Area, 

Nose  Gear 

11-20 

F10.3 

PINDN 

Net  Orifice  Area  at  STROKN,  Nose  Gear 

(in2) 

Card  15 

]-5 

15 

NSCM 

Number  of  Area  Changes  in  Main  Gear  Metering 
Tube 

*Card  16A,  16B, - 


1-10 

F10.3 

STR0KM 

Stroke  Corresponding  to  Orifice  Area, 

Main  Gear 

11-20 

**Card  18 A, 

F10.3 

18B . 

P1NDM 

Net  Orifice  Area  at  STROKM,  Main  Gear 

(in2) 

1-10 

F10.3 

SIMAIN1 

Mode  Shape  Deflection  Rear  Main  Gear 

11-20 

F1Q.  3 

SIMAIN2 

Mode  Shape  Deflection  Front  Main  Gear 

21-30 

F10.3 

SIN0SE 

Mode  Shape  Deflection  Nose 

31-40 

F10,3 

SICG 

Mode  Shape  Deflection  CG 

41-50 

F10.3 

SITAIL 

Mode  Shape  Deflection  Tail  Sfafion 

51-60 

F10, 3 

SIPS 

Mode  Shape  Deflection  Pilot  Station 

*  One  card 
program 

is  required  for  each 

stroke-metering  pin  combination  read  into  the 

**  One  card 

is  required  for  each 

flexible  mode. 
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Runway  Profile  Magnetic  Tape 

The  runway  profile  is  read  into  the  program  from  a  magnetic  tape, 
format  for  this  tape  is  shown  in  Table  V. 


Card 

Column 

Format 

RUNWAY 

Variable 

Name 

TABLE  V 

PROFILE  MAGNETIC  TAPE 

Definition 

Card  1 

• 

1-80 

8A10 

SITE 

Runway  Profile  and  Direction 

Card  2 

1-6 

16 

NPTSS 

Number  of  Runway  Elevation  Points 

*Card  3,  A, 

•  •••«»  Ni*2 

1-70 

10F7.  3 

ELEV 

Runway  Profile  Data 

*  One  card  required  for  every  ten  runway  profile  elevation  .-'oints. 


The 
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TABLE  VI 


Sample  Aircraft  Data  Deck,  KC-135 


XC-135  273000  POUNDS 


273000. 

46.0 

502.0 

57772000. 

581.0 

500.0 

1.00 

55000. 

278.88 

0.603 

2433.0 

0.06 

2634.0 

343.0 

2.0 

1.0 

1317.0 

24.51 

64.77 

82.3 

150.0 

222.0 

392.2 

1810.0 

3.149 

92.0 

92.0 

42352.0 

11957.0 

0.01 

0 

5 

0.001 

1.064 

6.41 

1.129 

9.86 

1.129 

13.88 

1.160 

16.58 

1.240 

3  3 

0.10 

1.483 

12.13 

1.483 

22.00 

1.980 

8 

4.5 

7.0 

4.5 

10.0 

8.0 

0.0 

1.5 

0.0 

1.0 

1.0 

0.0 

4.0 

0.0 

2.5 

4.0 

-3.0 

-3.0 

-3.0 

2.5 

-3.0 

24.2 

-68.0 

23.0 

-.410 

-88.0 

7.0 

-5.5 

7.0 

0.0 

-10.0 

-1.0 

-16.0 

-1.0 

-21.8 

-24.0 

-5.0 

15.0 

-6.0 

3.5 

30.0 

81903.4 

10.62 

946.62 

14.51 

7i.28.17 

15.58 

34028.44 

18.47 

724890.1 

20.80 

48741.8 

23.69 

121052.2 

29.47 

63481.19 

36.88 

1.227 
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Deck  Setup 

Figures  1,  2,  3,  and  4  contain  schematic  diagrams  of  the  source  deck 
setup  for  conventional,  F-4,  F-lll,  and  C-5A  aircraft  simulations. 

For  conventional  aircraft,  no  modifications  of  the  basic  deck  are  re¬ 
quired.  Only  the  correct  aircraft  input  data  must  be  used.  This  basic  source 
deck  setup  is  shown  in  Figure  1. 

In  order  to  simulate  the  F-4,  however,  changes  must  be  made  to  the  basic 
program.  This  change  is  a  result  of  the  main  landing  gear  having  a  double¬ 
acting  strut,  which  is  described  in  Volume  1  of  this  report.  The  F-4  source 
deck  is  formed  by  removing  the  subroutines  Taylor  and  IC  from  the  basic  source 
deck  and  replacing  them  with  the  Taylor  and  IC  subroutines  designated  for  use 
with  the  F-4.  Thi3  is  shown  in  Figure  2.  The  modifications  of  the  aircraft 
data  cards  described  in  this  volume  must  also  be  made. 

Changes  to  the  basic  source  deck  for  an  F-ill  simulation  are  shown  in 
Figure  3.  The  use  of  an  articulated  main  landing  gear  and  a  fluted  metering 
pin  on  this  aircraft  necessitates  this  change.  The  F-lll  source  deck  is  formed 
b>  removing  the  Taylor  and  IC  sun routines  from  the  tasic  source  deck  and  replacing 
them  with  the  Taylor  and  IC  subroutines  designated  for  the  F-lll  aircraft.  The 
aircraft  input  data  changes  shown  in  this  volume  must  also  be  made. 

Due  to  the  complexities  of  adding  an  extra  set  of  main  landing  gear  to 
the  basic  computer  code  by  substitution  of  subroutines,  a  completly  new  deck  is 
provided  for  the  C-5A  and  other  multiple  strut  aircraft.  The  source  deck  setup 
for  the  C-5A  aircraft  is  shown  m  Figure  4.  The  modification  of  the  aircraft 

data  input  for  the  C-5A  is  shown  in  this  volume. 
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\T^ 


SUBROUTINE  IC 


/  SUBROUTINE 
TAYLOR 


/67H9 


/AIRCRAFT  DATA 


/SUBROUTINE  TLOOK 


SUBROUTINE  COEF 


/ - 

/SUBROUTINE  IC 

/  JLcj-a) 

/  /  / 

'  X SUBROUTINE  1  , 

''PROGRAM  TAXI 


^/LGO 

^/lAP(ON) 

'SETCORE. 


^/NOREDUCF.. 

^/request 

""JOB  CARD  I 


FIGURE  2.  Source  Deck  Setup  For  F-4  Aircraft 
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SECTION  III 


PROGRAM  DESCRIPTION 

Program  Flow  Charts 

The  following  pages  contain  flow  charts  of  the  program  TAXI.  The  basic 
computer  program  is  flow  charted  entirely.  The  program  is  broken  down  into  its 
individual  routines  and  each  routine  is  flow  charted  separately.  The  flow 
charting  symbols  and  their  definitions  are  shown  in  Figure  5.  The  conventional 
direction  of  flow,  from  top  to  bottom  of  the  page  is  used. 
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PROGRAM  FLOW  CHART  SYMBOLS 


SYMBOL 


DEFINITION 


Operation  Box 


Card  Input 


Tape  Input 


Printed  Output 


Decision 


Subprogram  Execution 


Program  Statement  Number 


Page  Connector 


CZD 


Termination 


FIGURE  5.  Flow  Chart  Symbols 


TAXI 


TAXI  is  the  main  routine  which  directs  the  entire  program  in  the  sequence 
of  operations  and  the  calculations  to  be  made.  It  reads  aircraft  data  from 
cards  and  the  runway  profile  from  a  magnetic  tape  and  outputs  this  dara.  The 
runway  profile  is  both  normalized  to  the  first  elevation  point  and  detrended  in 
TAXI.  TAXI  calls  the  subroutines  IC,  COEFF,  and  TAYLOR.  IC  returns  aircraft 
initial  conditions  which  are  used  to  initialize  aircraft  dynamic  variables. 
COEFF  returns  the  coefficients  of  the  polynomial  fits  to  the  runway  profile 
segments.  Taylor  returns  the  solution  to  the  differential  equations  of  motion. 
TAXI  then  determines  if  this  data  is  to  be  printed  and/or  stored  for  use  in 
the  Calcomp  plot.  TAXI  also  directs  the  Calcomp  plotting. 
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WRITE  OUTPUT 
COLUMN 
HEADINGS 
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26 


INITIALIZE 
2ND  DERIVATIVES 
OK  FLEXIBILITY 
GENERALIZED 
COORDINATES 


COMPUTE  TOTAL 
ACCELERATION 
DUE  TO 
FLEXIBLE 
MODES 


COMPUTE  TOTAL 
ACCELERATION'S 
SUM  OF 
FLEXIBLE 
AND  RIGID 
BODY  MOTIONS 


HORIZONTAL 
VELOCITY 
GREATER  THAN 
\  TAKEOFF  y 


f  AIRCRAFT  \ 
EXCEEDED  RUNWAY 
V  SEGMENT  / 
\  LENGTH  / 


es 


( 


Taylor  is  the  subroutine  which  computes  the  dynamic  motion  of  the  aircraft. 
Using  runway  profile  information  and  aircraft  data  from  TAXI  and  net  orifice 
areas  by  calling  subroutine  TLOOK  it  calculates  forces  in  the  landing  gear  and 
tires.  Subsequently,  the  aerodynamic  and  mechanical  forces  and  moments  acting 


on  the  aircraft  fuselage  are  found.  Taylor  then  sets  up  the  differential 
equations  of  motion  and  integrates  them  using  a  three  term  Taylor  series  method. 
The  dynamic  variables  are  returned  to  TAXI  for  printing  and  storage. 
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SUBROUTINE  TAYLOR 


1 


o 


COMPUTE 
AERODYNAMIC 
DRAG  &  TIRE 
DRAG 


YES 


. SET  THRUST 
EQUAL  TO 
TOTAL  DRAG 
FOR  CONSTANT 


SPEED 
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Subroutine  IC 


This  subroutine  calculates  the  static  conditions  of  the  aircraft  needed  to 
initiate  an  aircraft  simulation.  From  statics,  the  reaction  force  at  the  main 
and  nose  jeats  are  calculated  by  summing  moments  and  forces.  Tire  deflections 
and  main  and  nose  gear  strut  stroke  are  computed  from  these  reactions.  Using 
these  values,  the  initial  aircraft  pitch  angle  and  the  initial  vertical 
position  of  the  CG  are  found.  These  values  are  returned  to  the  main  routine 
TAXI. 
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SUBROUTINE  IC 
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Subroutine  COEFF 


COEFF  calculates  a  third  order  polynomial  fit  to  a  four  foot  runway  pro¬ 
file  segment.  Given  three  runway  elevation  points  and  the  slope  from  the  end 
of  the  previ-us  runway  segment,  a  set  of  four  simultaneous  equations  is  solved. 
This  solution  yields  the  four  constant  coefficients  of  the  runway  segment 
polynomial.  This  subroutine  calculates  a  runway  profile  segment  fit  for  each 
s-.t  of  landing  gear  after  each  four  foot  traversal  of  the  runway.  The  coeffi¬ 
cients  are  returned  to  TAXI  and  used  in  Taylor  for  computing  the  runway  profile 
elevations  at  each  time  step  during  a  simulation. 


SUBROUTINE  COEFP 


COMPUTE 
COEFFICIENTS 
OF  3rd  ORDER 
POLYNOMIAL 
CURVE  TO  FIT 


THROUGH 

RUNWAY 
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Subroutine  TLOOK 


TLOOK  performs  a  linear  interpolation  of  the  values  in  the  table  of 
metering  pin  diameter  vs.  stroke.  For  an  aircraft  witn  a  metering  tube  or 
fluteu  r.c.tering  pin  ,  this  table  becomes  net  orifice  area  va.  stroke.  In 
either  case  TLOOK  is  called  by  Taylor  and  furnished  with  a  landing  gear  strut 
stroke.  It  does  a  table  look-up  and  fits  a  straight  line  between  the  two 
points  which  straddle  thegiven  strut  stroke.  The  slope  and  y-intercept  of 
this  line  are  returned  to  Taylor  where  a  net  orifice  area  is  computed.  TLOOK 
is  called  at  each  time  step  during  a  simulation  for  eac-  i'et  of  landing  gear 


of  the  aircraft. 


SUBROUTINE  TLOOK 


RETURN 


SECTION  IV 


PROGRAM  OUTPUT 

The  results  oi  Ll»c  program  are  presented  in  two  forms.  These  are  a 
listing  of  ten  dynamic  aircraft  parameters  and  a  Calcomp  plot. 

The  listing  of  the  dynamic  aircraft  parameters  occurs  at  specified  time 
intervals  during  an  aircraft  simulation  which  are  larger  than  the  integration 
step  size.  For  a  take-off  simulation  printing  occurs  at  .01  second  intervals. 

If  a  constant  speed  taxi  simulation  is  made,  printing  occurs  every  .02  seconds. 
These  intervals  are  model  simulation  times  not  real  time.  Thus,  using  an 
integration  step  size  of  .001  seconds,  ten  or  twenty  solution  integrations  are 
made  between  every  line  of  printing.  The  ten  aircraft  parameters  selected  for 
printing  are  main  gear  stroke,  nose  gear  stroke,  main  gear  force,  nose  gear 
force,  speed  of  the  aircraft,  distance  down  the  runway,  tail  acceleration, CG 
acceleration,  pilot  station  acceleration,  and  simulation  time.  Other  variables 
may  be  printed  out  by  modifying  the  WRITE  statement  in  the  main  routine  TAXI 
and  putting  the  variable  in  COMMON' between  TAXI  and  the  routine  in  which  it  is 
defined.  A  sample  pe»e  of  printed  output  is  shown  in  Figure  6. 

The  other  form  .  ,•  put  is  the  Calcomp  plot.  On  this  plot,  CG  acceleration 

and  pilot  station  acceleration  time  histories  are  displayed  along  with  the 
runway  profile  time  history  as  seen  by  the  nose  gear  of  the  aircraft.  Aircraft 
speed  and  distance  and  runway  markers  are  also  plotted  at  specified  intervals. 

A  more  complete  description  of  the  Calcomp  plot  is  contained  in  Volume  I  of 
this  report.  A  photographic  reduction  of  a  typical  Calcomp  plot  is  shown  in 
Figure  7. 
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FIGURE  6.  Program  Output  Listing 
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FIGURE  7.  Typical  Calcomp  Plot 


APPENDIX  I 


PROGRAM  LISTINGS 

A  complete  listing  of  the  program  source  deck  for  conventional 
landing  gear  aircraft  and  the  C-5A  aircraft  is  given  belov.  Listing 
of  the  TAYLOR  and  IC  subroutines  for  the  F-4  and  F-lll  simulation  are 
also  included.  The  programs  are  separated  into  individual  routines 
for  clarity. 
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PROGRAM  TAXI (INPUT  ♦0UTPUT,TAP£5=INPUT,TAPE6=0UTPUT,TAPE3,TAPc2, 
1  TAPE?) 


♦  #*«*♦**♦***♦****•**♦*#♦**♦*#*■*•#**, ,*•»**♦*#***  **♦***♦*,>*  *#**♦*  ********** 

C  THIS  PROGRAF  MILL  PERFORM  M  SYMMETRIC 
C  OYNAMC  TAXI  ANALYSIS  ON  A  FLEXIBLE 
C  VEHICLE  TRAVERSING  A  GIVEN  RIGIO  RUhNAY 
C  PROFILE. 

**«, «*«,*«*«* o ****************** *9 ************ ******** *«***«*t  *« ******** 


L0Mt.eN/FLcXi/SIMAINU5)  *SIhOSE(15)  ,SICG(  IS)  ,SITAIL(15),SIPS(15) 
C0MM0N/FL£X2/NFM,GH(15) ,CMEGA(15) 

C0HMCN/FL£X3/Q(i5),Q0 (15) ,000(15) 
C0MM0N/X1/M,MM,MN,HCG,MM,HN,A,8,MMI 
CUMHCN/XZ/PACM, VCM»AAM, AHM.OAM 
C0MM0N/X3/PACN, VON, A AN« AhN.OAN 
COMMON/ X<»/SXM,SXN,SLH,SLh,TSN  ,TSN 
COMMGN/X 5/CL  «C0« AREA «THRL ST 
COMMON/ Xft/Z,REACTM,R£ACTh*CX,NTRUN,SP£EO 
COMMON/ X7/AM,8M, CM,OP,AN, BN,CN, Oh 

CQMMCN/X6/STROKM(20) ,PINCM(2Q) ,STROKN(20) ,PINON(20) ,NSCM,NSCN 
C0MM0N/X9/FSh,FSh,FTh,FTh,XMAlN,XK0S£, VELM,V£Lh 
COMMCN/X 1 E/7PM, 2PN 
DIMENSION  PLANE (6) , SITE ( fi) 

DIMENSION  T  1 12) *  T0(12t,  T00(12) 

DIMENSION  £L£V(7  GOO) 

DIMENSION  YP (**) , YPN(N) 

CIMENSION  CG ACC (160  3)  *  TIMEU830),  PSACC(iSOO),  TIMEl(lfcOO) 

C IMExSICN  PRCFU80C)  ,RMARK(20) 

DIMENSION  SSPLOT  <3Q3) ,  S T IME ( 3 0 C ) , OOPLOT ( 30 C ) 

REAL  KCG,KN,NN,MHI 
CALL  PLGTS(1G0«,QUH,3) 

PSA=u.O 

CGOUT=O.C 

rOR=0. 

x=a. 

LLL  =  1 
LL  =  0 
STOREl  * 

STCRE2  =  C. 

ST0RH3  =  0. 

STORE*,  =  C. 
hN  =  0 
II  =  0 
FP  =  1.0 
M  -  l 
ITT  =  G 


C  READ  ANO  PRINT  INPUT  CATA 


C  M* VEHICLE  HEIGHT  AT  CG  (POtNCS) 

C  A=C ISTANCE  MAIN  GEAR  TO  CG  (INCHES) 

C  9=C ISTANCE  NOSE  GEAR  TO  CG  (INCHES) 

C  MNI=NASS  MOMENT  OF  INERTIA  (L8  IN  SEC  SQ) 

C  PLAhE=  AIRPLANE  BEING  SIMULATED  ANO  GROSS  MEIGHT 
C  PSARM  -  UISTmNCE  OF  PILOT  STATION  TC  CG 
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C  TAILRH  =  DISTANCE  OF  TAIL  STATION  TO  C6 
C  TAKCFF*  TAKE-OFF  SPEED  (FEET/SEC) 

C  SPEEC  ^INITIAL  VEL  OF  AIRPLANE 
C  THURST=  TOTAL  AIRPLANE  THRUST 
C  CL-LIFT  COEFF. 

C  AREAsNlNG  AREA 
C  CD3CRAG  COEFF. 

C  UN-LEIGHT  OF  MAIN  GEAR  (EACH) 

C  NN=KEIGHT  OF  NOSE  GEAR 
C  SXH=  NUMBER  CF  PAIN  GEAR  STRUTS 
C  SXN=  NUH8ER  OF  NOSE  GEAR  STRUTS 
C  AH**  HV0RAULIC  PISTON  AREA  NCSt  SQ  INCHES 
C  AAN  PNEUMATIC  P1STCN  AREA  NCSE  SQ  INCHES 
C  AHH  HVORAULIC  PISTON  AREA  MAIN  SQ  INCHES 
C  AAM  PNEUMATIC  PISTON  AREA  MAIN  SQ  INCHES 
C  PAON  NOSE  STRUT  PRELCAC  PRESSURE  PSI 
C  PAOH  MAIN  STRUT  PRELCAC  PRESSURE  PSI 
C  VON  NOSE  STRUT  INITIAL  VOLUME  CU.  IN. 

C  VON  MAIN  STRUT  INITIAL  VOLUME  CU.  IN. 

C  OA.i  ORIFACE  AREA  MAIN 

C  CAN  ORIFACE  AREA  NOSE 

C  SLN=MAIN  GEAR  STRUT  LENGTH  UNLOADED  INCHES 
C  DISTANCE  FROM  CL  OF  AXLE  TO  CG  LINE 
C  SLN*NOSE  GEAR  STRUT  LENGTH  UNLOADED  INCHES 
C  DISTANCE  FROM  CL  OF  AXLE  TC  CG  LINE 
C  TSM  MAIN  TIRE  SPRING  CONSTANT  PER  STRUT 

C  TSN  NOSE  TIRE  SPRING  CONSTANT  PER  STRUT 

C  9X*TIHE  STEF  SIZE 

C  REAC  METERING  PIN  DESCRIPTION  STARTING  AT  ZERO  STROKE 
C  NSCN*5  OF  METERING  PIN  CHANGES  NOSE  GEAR 
C  NSCH==  OF  METERING  PIN  CHANGES  MAIN  GEAR 
C  NFM  =  NUMBER  OF  FLEXIBLE  NOSES 
C  SIXXXXC I)  x  HOOE  SHAPE  OEFLECTION  (NON  OXM.) 

C  SH(I)  =  GENERALIZED  MASS  CPCUNOS-SEC  SQ/IN) 

C  OMEGA  (I)  «  MODAL  FREQUENCIES  (RAO/SEC) 

READ  (5.1)  PLANE 
1  FORHATUAIO) 

REA0C5.5)  M»  A.B.MHI 
5  FORMAT ( 3F 10. l.FIZ.b) 

READ (5.10)  PSARM.TAILRM 
10  FORMAT (2F10. 2) 

R£A0(5»15)  I  PEEO.THRLST.TAKOFF 
IS  FORMAT(3F10.3) 

READ  (5.20  CL, AREA, CO 
20  FORMAT C3F10. A) 

REA0(5*25)  MM.HNf SXM.SXN 
25  FORMAT(Af 10.2) 

R£AO(St30)  AHN.AAN.AHMt AAM 
30  FORMAT  (•♦»  10.5) 

READ (5. 35)  PAOM.PAON, VON ,VCH .0AM.0AN 
35  FORMAT (6F10. 5) 

fi£A0(5«NG)  SLM.SLN 
NO  FORMAT (2F10.3) 

REAu(S.Ng)  TSM  .TSN 
*»S  FORMAT  (2F10.1) 

R£A9(5f 50)OX 
REAO (5.51) IFPLOT 
50  FORMAT (F10.3) 
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fi£A0(5,51)  NSCN 

51  FORMAT (15) 

fi£AQ(5»52)  (STROKN(I),PINOMI) « 1=1, NSCN) 

52  F0RMAT(2(F1<)*3>) 

R£A0(S*51)NSCH 

RcAO  (5*52)  (STROKH(I)  ,PINOP(I)»I=l,NSCM) 
fi£A0(5,51)  NFM 

R£A0(5»53)  ( SIHAIN ( I ) , SIKOSE (I ) ,SICG(I),SITAIL(I)*SIFS(I) 

1  * 1=1, NFM) 

53  FORMAT (5F10 • 3) 

R£A9(5,54)  (GM(I) ,  CK-GAd)  ,1  =  1,  NFM) 

54  F0R**AT(F15«2»F1Q.3I 
PM=WM/38€. 

MN=WN/386. 

MCG=W/  38  6 » 

TYPRUN  =  .01 
NTRUN  =  1 

IF (SP££0  .GT.  5.)  56,57 

56  TYPRUN  =  .01 
NTRUN  =  2 

57  WRITE (6, 55) 

55  FORMAT!  IK.  46X,  **♦♦•♦*♦*•*••****♦*♦♦*♦**•♦♦***♦*♦**•*****) 

WRIT£(6,6C) 

6C  FORMAT (46X,46H******##,M,#*4*#**  INPUT  OATA  **♦*#****♦**♦*♦*♦) 
MRITE(6,6l) 

61  FORMAT!/// IX, 41M*******4*  GENERAL  AIRCRAFT  OATA  ****♦*•*♦) 

WRITE (6, 62) PLANE 

62  FORMAT!//, IX, 8A10) 

WRITE (6 ,65)  W*MM«MN,A«B*PHI 

65  FORMAT !///,5X,*N=*,F10.1,5X, F10. 2, 5Xt*MN=*,F10. 2, 5X, 
1#A**,F10 .3,5X,*B=*,FiQ.3,5X,*HMI=*,FlE.fl) 

WRITE (6,66)  S XM , S X N, SLM, S LN , PS ARM , TA ILRM 

66  FCRMAT!//,5X,*SXM**,Ftt.l,5X,*SXN**,F9.1,5Xf*SLM=*,F9.1,5X, 
l*SLN=*,F8.1,5X,*PSARM=*,F7.1,4Xt*TAILRM=*,F7.1) 

WRITE (6,70)  AAM, AHM,PAOM ,VCM,OAM,TSH 
70  FORMAT !//,5X »*AAP=* ,F8»2 ,5X**AHH=*»F8. 2, 6X» *PAOH=*,F9. 2»4X» 

1* WOM=*,F8 .2, 5X»*CAM**,F8  «2,5X, *TSM  =*»F10.2) 

WRITE (6,74)  AAN, AHN, PAON, VON,OAN,TSN 

74  F0RHAT(//,5X,*AAN**,F8.2,5X,*AMN**,F6.2,6X,*FA0N=*,F9.2,4X, 
1*¥GN=*»F6.2,5X,*0AN=*,F8.2,5X,*TSN  =*,F10.2) 

WRITE (6, 75)  CL, CO, AREA, SPEEO, THRUST, TAKOFF 

75  FORMAT (// *5X »*CL=*  ,F9.3, 5 X,*CQ=* ,F9. 3,6X,*AREA*» ,F8 .2, 5X, 
1*SPEEOs*,F7.1,4X.*THRUST»*,F7.0,3X,*TAKOFF«*,F10.2) 

WRITE  (6,76)  (STROKMI)  ,PINON(I)  ,1=1, NSCN) 

76  FORMAT (//,3X,*STR0KE  NOSE  PIN  CIAMETER  *//,2 (3X,F10 .3) ) 

WRIT £(6,77)  (STROKPC I), PINDM(I), 1=1, NSCM) 

77  FORMAT!//, 3X,*STR0K£  MAIN  PIN  CIAMETER  *//,2(3X,F10.3) ) 

WRITE (6, 60) 

80  FORMAT (1P1,1X**PC0E  SIPS  SIN0SE*,4X,*  SICG  SIMAIN  SITAIL 

IOMEGA  GEN.  PASS*) 

CO  65  1*1, NFM 

65  WRITE (6,50  ) (I, SIPS! I) ,SINCSE(I) «SICG( I), SIMA IN (I), SITAIL! I) , 
IOMEGA (I) ,GM( I) ) 

90  FORMAT!//,  2X, I2,6(F7.2,2X) ,F10 ,1) 

C  ZMI  *  INITIAL  MAIN  GEAR  POSITION  INCHES 
WRITE(6,S5) 

95  FORMAT (///,46XvA0H*********  INTIAL  CONDITIONS 
C  ZNI  ^INITIAL  NOSE  GEAR  POSITION  INCHES 


C  ZCGI*INITIAL  C.G.  POSITION  INCHES 
C  THETAX*XNXTX«L  PITCH  ANGLE  DEGREES 

C  THE  ABOVE  PARAMETERS  ARE  CALCULATED  IN  SUBROUTINE  IC 
CALL  IC(  ZCGI*ZMI*ZNX*TH£TAX) 

NRITE (6* 10$)  ZHI.ZNItTHETAl.ZCGI 

105  FOAHAT(//.5X,6MZHIs,Fi0.3*5X.6HZNIs,F10.3*5X*7HTHETAIs*F10.6* 
15X«5HZC6I>*F10*3) 

REACT N« ( ZNI*TSN  )*SXN 
REACTH»(2NI*TSH  )»SXH 
NRITE (6.106)  XHA IN* XNOSE  *  REACTN*  REACTH 

106  FORMAT  (/*3X«*XHAINs»« F0.3.3X**XNOSE«**FO*3*3X**REACTNs*,F10.0, 
13X»*REACTM** *F10  *0) 

C  READ  RUN NAY  PROFILE  OATA  (ELEV) 

C  SXTEs  RUNHAV  PROFILE  ANO  DIRECTION 
C  NPTSSs=  OF  RUNNAY  ELEVATION  OATA  POINTS 
NRITE  (6*111) 

111  FORMAT (*T*I 
00  112  I«l*50 

112  ELEVdlsO.O 
READ (2*1)  SITE 
REA0(2«116)  NPTSS 
NPTSS*NPTSS*50 

116  FORHAT(IS) 

L0«51 

117  LSD  s  LO  ♦  9 

REA0(2* 116)  (ELEV(I) *I*L0*LS0) 

110  FORMAT (10F7«3) 

MO«LO 

IF (LS0*GE tNPTSS)  GO  TO  120 

LOsLDMO 

GO  TO  117 

120  £LEV1s£LEV<S1) 

00  125  Is51«N10 
125  ELEV(I)  s  ELEV(I)  -ELEV1 
CALL  RETURNS (2) 

OISTANsO. 

LSOsLSO-OO 
SLP*(LSO*50) *2 
SLPs£LEV(LSO)/SLP 
00  126  X«51  *N10 
ELEV(I) *£L£V (I) -SLP*CXSTAN 
126  0ISTAN«0I$TAN42* 

IVALs(A*C)/26. 

NRITE (6* 130) 

130  FORHATdNifSX**  RUNNAY  PROFILE  OATA  NQRMALIZCO  (SLOPE  REMOVIO)* 
150X«*F|£T  00 NN  THE  RUNNAY*) 

NRITC(Otl)  SITE 
L01«l 

139  LS01  ■  L01  ♦  9 
LPRIN  ■  LSOi  •  t 

NRITE (6*160)  (ELEV(I) *S*L01*LS01) * LPRIN 
160  FORMAT (IX «10F10*S*10X«IO) 

L01SL01*10 

IF(L01*GE«N10)  €0  TO  169 
GO  TO  139 

169  NRITE (6* 196)  L01 

190  FORMAT (62N******9***  CNO  RUNNAY  OATA  L01«*I9) 


51 


I»R  IT  £  (7*1)  SITE 
ENORUN  =  NPTSS*2 
FRK= (ENORUN- 100 • )/10C0 * 

RM=(((ENORUN-100.) -FLOAT  (MRMMOOO)  )/2.  DUOQ .  - 
EH0RUN=£K0RUA-90. 

T (2)  *2CGI' 

T0(2)  *  0. 

T(4)  =  ZMI 
T0(4>  a  C. 

TC6»  =  ZM 

T0(6)  a  c.  *  ».  ,  .•••.«. 

T ( 41  =  TNETAI 

lo*8!  =  C;  ...... 

1C18I  a  C. 

TO(10)  a  SPEEO 
NRIT£(6» 155) 

155  FORMAT ( 1H1#  7X**  XMAIN  XNOSE  FSM  FSN 

1  TAILAC  SPEEO  OIST.  CGACC  PSA  TIM 

2E*) 

cz=o. 

ZPMaO. 

ZPN=Q. 

ZDOT=Q. 

ZNOOT=0 • 

J*1 

160  YP(1)  =  ZOOT 

YPC2) *ELE VC J) 

YPC3)aELEtf<J4l) 

YP(A)*ELEVtJ42) 

YPN(l)aZAOOT 
VPN(2)s£LEV( IVAL) 

VPN(3)*£LcV(IVAL«l> 

YPNU)*ELEV(XVAL*2) 

ZaOZ 

CALL  GOEFF  ( YPN* AN*3N,CN*0N) 

CALL  COEFF  (YP*AN,BH,CH*CH) 

ZOOTaBH 

ZNOOYaBN 

IF CH.GT • 1)  GC  TO  245 
M  =  M  ♦  1 

165  WRITE (6*170)  XMAIN*XNOS£«FSM*FSN«TAILACtTO(iO) *T(lCt * 

1CGOUT  ,PSA*X 

NRIT£(7*  950)  TAILAC »CGOUT* PSA 
900  FORMAT (3X*F10»3) 

TIMEX  a  C. 

170  FCRMAT(3X»2(3X*F10.3)«2I3X*F10.0)«5(3X*F10.3) *2X*F6.2) 


IF(IFPLOT.EQ.l)  GO  TO  230 

IF  (ABS  (CGOUT  )  .LE .A9S  (STCRE2)  . ANC. ASS f ST0RE2)  .GE.ABS (STORED  ) 
ICO  TO  180 

175  STORE1  =  ST0RE2 
ST0RE2  a  CGOUT 
GO  TC  19 C 

180  IF ( X-TIME  (NN) *GT • *08)  GO  TO  185 

IF (A9S(ST0RE2) .GT..300CC)  GO  TO  185 
GO  TO  175 
165  NN  MM  1 

CGACC (NN)  >  STORES 
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PROF (NN)  =  ZPN 
TIPF.(NN)  *  X  -  10.  *  OX 
IF(TIHE(KN).tT.O.)  TINEJRN)  =  . 

STORE1  =  ST0RE2 
ST0RE2  =  CGOllT 

190  IF  (A6S  IPSA). LE.A6S (STORED) .ANO.ABS (STORE 4) .GE.A6S (ST0RE3)) 

1G0  TO  205 

200  ST0RE3  -  ST0RE4 
STORE4  =  PSA 
GO  TC  215 

205  If  (X-TIfcEiaL).GT..06)  GC  TO  210 
IF (A6S(ST0R£4) .GT« .32)  GC  TO  210 
GO  73  200 
210  LL  =  LL  ♦  1 

PSACCtLL)  =  ST0RE4 

TIF£1(LL)=X  -  10*0X 

ZF(TIHEKLL)  .LE.O.)  TIMEl(LL)  =  .001 

ST0RE3  =  ST0RE4 

STORE1*  =  PSA 

215  IF (A8S(T (10) -RN) . LT.5.C)  GC  TO  220 
GO  TO  225 
220  II  *  II  ♦  1 

RHARK(II)  =  X 
RH  =  RH  ♦  1000. 

ITT  =  II 

225  IFU8SU-PP)  .GT.  .005)  GO  TO  23G 
SSPLOT(LLL)  =  TO (10) 

STIHE(LLL)  =  X 
COPLOT (LLL)  =  T(10) 

LLL  =  LLL  ♦  1 
PP  *  PP  ♦  lit 

n(i«flfl««*«f«««*«««»4H«**««*M*««««*»*M*«*f***»*«***********«*********** 

230  IF (HOR.LE .00 •)  GO  TO  235 
HQR=0 . 

4RITE(6»155) 

235  FOR  =  HDR  ♦  1. 

240  X  =  X  *  CX 

CALL  TAYLOR CT « TO. TOO) 

IF  CT (10) .GE. E NORUN)  GO  TC  260 
IF (X.GE.300. )  GO  TO  263 
TIPEX=TIPEX*GX 
QOOCG=C • 

G00PS=0. 

GOOTAL*0  . 

00  241  I=1»NFM 

QDOTAL=COOTAL+QOO(I)*SITAIL(I) 

000CG=Q0CCG*000( I) *SICG(I) 

241  QOOPS*QCQPS*QOO(I)*SIPS(I) 

TAILAC=(T00(2)*TAILRP*T00(<J)  )/366.  ♦Q00TAL/366. 
PSA=(T0DC2)-PSARP#T00C8))/3«6.  ♦QQOPS/386. 

CG0UT=T0C(2)/336.  4Q00CG/366. 

IF (TO(IO) .GE.TAKOFF)  GO  TO  250 
IFIZ.LT.4.)  GO  TO  245 
C2=Z-4. 

J=J*2 

IVAL=l¥AL*2 
GO  TC  16C 

245  IF (TIMEX.LT.TYPRUN)  GO  TO  24G 
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GO  TC  165 
250  WRITE (6*  255) 

255  F0RMAT(5X,*  THE  VEHICLE  HAS  TAKEN  OFF*) 

260  WRITE(6,262)  ENDRUN* T (10 ) 

262  FORMAT <3*,*  ENO  OF  RUNNA Y*,2F10 . 3) 

TOTAL  =  X/TYPRUN 
WRITE (6, 3G0)  TOTAL 

300  FORMAT {  *  TOTAL  NUMBER  OF  POINTS  ON  TAPE  IS*  F12.0) 

263  MGH  =  X 

IF (IFPLOT .EO.l)  GO  TO  265 
XLONG  =  FLOAT(MGM) 
np:te(6,2*5)  NN,LL 
265  F0RMATC2I20) 

IF(NN.LE.180C.QR.LL.LE.1800)  GO  TO  266 
hRIT£(6,267) 

267  FORMAT (3Xt*TME  ARRAYS  CGACC  OR  PSACC  OR  PROF  HAVE  EXCEEOEO 

i  their  dimensioned  size*) 

266  CONTINLE 

C266  CALL  PLOTS(OATA,438) 

CALL  FACTOR  (2.0) 

CALL  PLOT(3«C,-li.e,-3) 

CALL  PLOT (3. « .7,-3) 

TIME (NN*1)  =  0.0 
T IME  (NN*2)  =  1.0 
T IME1 (LL+1)  =0.3 

TlHcKLL  +  2)  =  1.0 
CGACC (NN  +  1)  =  -1.0 
CGACC(NN*2)  *  1.0 
PSACC (LL*1)  =  -1.0 
PSACC (LL*2)  a  1.0 
CALL  SCAL£(PROF* 10  ««NN,1) 

PROFIO  =  PROFCNN+2)  *10. 

IF (PROF 10. 6T .10.5)  GO  TO  270 
PRCF (NN+2)  =  6.0 
IF (NPTSS.GE.lflOO)  GO  TO  270 
CO  269  1=1, NK 

269  PROF (I)=PROF (I) *36. 

270  CALL  AXIS (0 • ,0  *«11HTIH£  (SEC.) » -11* XLONG .0.0 .TIME (NN+1) » 
1TIME (NN*2) *0 ) 

I XLONG  =  XLONG 

CALL  PLOT  (XLONG, 1.1,3) 

CALL  PLOT  (;. ,1.1,2) 

CALL  SYM6GL  ( .1, .4, .105, 8HOISTANCE,0 ., 8) 

CALL  SYMBOL  ( .1, .2, .10 5, 5HSPEEO, C. ,5) 

OO  275  1=1, IXLONG 

CALL  NUHEER(STIHE (I) , .2, • 1C5,SSPL0T (I) ,0  .0, 4HFA.0) 

275  CALL  NUHEER(STIME (I) , .4, .105, OOPLOT ( I) ,0 .0, AHF5.G ) 

CALL  PLOT  (XLONG, 1.5,3) 

CALL  PLOT  (0 .,1.5,2) 

00  280  I=1,ITT 

280  CALL  SYMBOL (RHARK(I) ,-*l«.2A5«70 ,0 .0,-1) 

CALL  PLOT  (XLONG, 1.9,3) 

CALL  PLOT  (0., 1.9,2) 

CALL  SYMBOL (A., 9.,. 1A, PLANE, 0*, AO) 

CALL  SYMBOL (A. ,8 .».1A*SITE»0.,AC) 

XL0NG2=XL0NG/2. 

C  CALL  SYMBOL (XL0NG2 ,9.».1A»11HAFF0L  -FYS-,0.,11) 


CALL  SYHBOLC XL0NG2, 3 .»  ,14»10HMPAFe  OHlQ»0.»lO) 

CALL  AXlSCo!!-l!o,3HCGA, 3,2.0 .90 t,C6ACC(NN*l»*CGACCCNN^2lt-i» 
CALL  NUMBER CO. 4 ,-.4,  .Iu5«-.4,Q0.C»4HF4.1) 

CALL  NUNEERC6.5,  .4,  .105,  ♦.4,04. 0  »4HF4. 1) 

CALL  PLOTCC. ,-l.Q»-3) 

CALL  LINECTIP£,CGACC»NN,1»Q,G> 

CALL  PLOT  CXLONG, 3.1,31 
CALL  PLOT  CO • ,3. 1,2) 

CALL  AXIS  CO. ,-1.0, 3HPSA, 3,2.0,90 PS ACC CLL+1) , PS ACC CLL+2) ,*1) 
CALc  PLOT C 0 . , C * , 3) 

CALL  PLOT  CXLCNG,0.,2) 

CALL  PLOT  CXLONG,. 4, 3) 

CALL  PLOT  C“ .1, .4,2) 

CALL  NUHEERC0.5, .4, .105, +.4, 00 . G.4HF4. 1) 

CALL  NUMBER  CO .4,“ .5, .105 ,-.4, 00 .u ,4HF4 .1) 

CALL  PLOTC0.,-1.0,-3) 

CALL  LINE CTIPtl,PSACC,LL ,1*0, 0 ) 

CALL  AXIS C-i.5^0. »15H£LE VAT ION  CIN.) ,15,10. 0,90.0, PROFCNNU) , 
lPR0FCNN+2),-l) 

CALL  PLOT  CO. ,Q.,-3) 

XPR0F=.25*ASSCPR0FCNN4l) l/PROF CNN*2> 

CALL  SYHEOL  C2.,XPR0F,.14,15HN0S£  GEAR  TRACK, 0., 15) 

CALL  LINECTIPE, PROF, NN,1, 0,0) 

XSTOP*XLCNG45. 

CALL  PLOT  CXST0P,Q.,-3| 

WRITE C6, 290) 

CALL  PLOTCO. ,0.,40) 

CALL  PLOTE 
STOP 

FORMAT C*S*) 

ENO 


SUBROUTINE  TAYLOR CT *TD* TOO) 

CCN.MCN/FLEX1/SIHAINC15) »SIKOS£C15) *3ICGC15) .SITAILC15) ♦ SIPS (15) 
COMMON/ F LcX2 / NFM *  GH (15) * CMcGA ( 15 ) 

C0MM0N/FIEX3/QC15>*QDC15>, 000(15) 

C0MM0N/X1/M*  MH*MN  *MCG*HM  *HN*A*  B*MHI 
C0MM0N/X2/PACM* VON* AAM* AHH*OAM 
C0MMCN/X3/PACN* VON* AAN* ANN«OAN 
COMMON/ X4/SXN*SXN*SLM*SLN*  TSM  ,TSN 
C0MM0N/X5/CL *C0«  AREA* THRUST 
COMMON/ X6/Z*R£ACTM* REACTS* 0X*N?RUN*SPE£O 
C0MM0N/X7/AN,Brt,CM*0M,AN*8N»CN,0N 

CUMMuW/ XS/STROMM  €201  ,PINCMC2Q> , STR0KNC20 )«PIN0NC2C) *NSCH»NSCK 

COMMON/ X9/FSM«FSS» FT N*FTS*XM AIN  *  XNOSE* VELH*  VELN 

COMMON/ XlQ/ZFMtZPN 

DIMENSION  T ( 12) *  T0C12),  T00(12> 

REAL  MCG*HM«MN«MMX 
1  Z=Z+T0C1C)#0X*T00 ClQ)*0X**2/2. 

12C  ZPH=AM*BM*Z*CM*Z**2+CM*Z**3 
ZPN=AN-»BK*Z*CN*Z**2*0N*Z**3 
QTN=0 • 

CTM*0. 

CTON=Q. 

QTOM*0. 

00  130  1=1* KFM 
QTN=QTN4Q(I )*SINCS£  CI) 

OTf=aTM*CtI)*SIMAIKCI) 

QTOM*QTCM*QCCI)*SIMAlNCI> 

130  QTON*GTON+QOCI) *SINOSEC 1) 

XNOSE  *  (TC2)  -  8  *  TC8)  -  TC6))  ♦QTN 
XMAIN  *  CT (2 )  ♦  A  *  TC8)  -  TU))  ♦QTM 
VELN  *  TC (2)  ♦  A  *  TOCS)  -  T0C4)  ♦QTOH 
VELN  =  TCC2)  -  B  *  TOCS)  -  T0C6)  ♦QYDN 
IFCXMAlN.GE.fi.)  XMAIN*-. 1 
IF CXNOSE.GE. C.)  XNOSE*-.  1 
IFC VELM.EQ.O .)  V£LM=-.l 
IFCVELN.EQ.C.)  VELN=-.l 
XHLK=ABS CXMAIN) 

XNLK=AQS  CXNOSE) 

C  NOSE  ANC  MAIN  OAMPING  COEFF 

CALL  TLOCKlXMLK.SLOPEM* VCEPM»STROKM»PINr  k » NSCM) 

CALL  TLOCKCXKLK«SLOPEN»YC£FN*STRCKN*PXN  ..  -  SCN) 

AOM  =  CAM  -C  CSLOFEM'XMLK+YCEPN)*^)*.//.  j9 
AON  =  OAK  - C  CSLC  FEN*XNLKt  YCEPN) **2) *• /S539 
CON=  C. COCOS'*  CA«N**3.))/C2.M.9*A0N)**2) 

COM=C.OOOC8* CAHM**3.))/C2.*C.9*ACM)**2) 

C  NOSE  ANO  MAIN  STRUT  PNEUMATIC  FORCES 
SSMMPAOMVOM)/  C  C  CVOM/AAM)  -XHLK)  ) 

SSN=  CI'AO'  *VOK)/  C  C  CVON/AAN)  -XNLX)  ) 

FTM  *  SXM  ♦  TSM  *  CTCA)  -  ZPM) 

FTN  =SXN*  TSN  *  CTC6)  -  ZPN) 

IFCFTM.GT.O.)FTM»Oo 
IFCFTN.GT .0 • )FTN=0 . 

FSMSXN*  C-SSK*C0N*VELN*A6SCVELN)  ) 

FSM=SXH*  C-SSM  ♦C0M*V6t.M*ABSC  VELM)  ) 

VLIFT  *• G31189*CL*Aft£A*  tTOf 10) #TC  CIO ) ) 

0RAGA*VLIFT*CD/CI 

CRAGT  =AeSC.C25*FTM*.025*FTN) 

IF  CNTRUN  . tQ. 1)  GO  TO  125 
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IF (TO <10  )  .LT.  SPEEO)  60  TO  125 
THRUSTs0fiAGA*0RA6T 

TOO (2)  =  («FSN-FSH-MCG*3e6.4VlIFT)/MCG 
T00(4)  =  <FSt'-FTH-386.*SXM*HN>/lHH«SXM> 
T00(6I=(FSN-F?N-HN*336.*SXNI/CHN*SXN) 

T00(8)  *• (FSF*A  -FSN*B  -ORAGT  * (SLM+XHAIN) ) /MHI 
TDD(IO)  -  (TFRUST-ORAGA-CRAGT) / C tMCG*12» > ) 

DQOOa»  =  -<SI^AlNCI)«<FSH-REACTH)*SINOS£(I)4  (FSN-REACTN) 
14.10*OH£GA(I l*QO (II^GHdUOHEGAf  II**Z*GH(I)*Q(I> )/GH(l> 
CO  1001  1  =  2,10,2 

▼ (I»  s  7(1)  ♦  TO(I)*OX  ♦  (TOO (I) #0X**2)/2» 

L  TO(i)  =  TO(I>  *  TOO(I)*OX 
00  1002  1=1, AFH 

Q(I)=Q(I)*QO(I>*OX4(QOOm*OX**2)/2. 

2  CO (I) XQ0 (I) +C00( I)*OX 
RETURN 
ENO 
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SUBROUTINE  IC(  ZCGI, 2HI.ZNI* THETAI) 

COMMON/ FLEXi/ SIM  A  INC  15)  *SlNOS£(15)»SICG(15)  ,SITAIU15)  ,  SIPS  (15) 
CCHMCN/FLEX2/NFH , GH ( 15 ) , OMEGA ( 15 ) 

C0MM0N/FLEX3/Q (15) »Q0(15) ,GOO < 15) 

COMMON/ Xl/M'NM«MN,MCGtMH*MN»A, 8 tMMI 
COMMON/ X2/PA0H* VOHf AAMf AHM«CAM 
C0MM0N/X3/PACN« VON«AAN» AHNtOAN 
COMMCN/X<t/SXK*SXN» SLMtSLN* TSM  ,TSN 
COMMCN/XS/FSMfFSNyFTMyFTNtXMAIN*  XNOSE*  VELHy VEEN 

C  THIS  PROGRAM  MILL  FIND  THE  INITIAL  CONDITIONS  FOR  TAXI 
C  rOR  mXH  MAI>  AND  SXN  NOSE  GEAR 
C  ZMI*  MAIN  GEAR  TIRE  OEFLECT ION  35  PERCENT 
C  ZNI=  NOSE  GEAR  TIRE  DEFLECTION  35  PERCENT  OEFLECTION 
C  ZCGI=  CG  OEFLECTION 
C  THETAI=  PITCH  ANGLE 

C  XMA IN=  MAIN  LANDING  GEAR  STATIC  STROKE 
C  XNOS£=  NOSE  GEAR  STATIC  STROKE 


RM=M/ (1* ♦A/BI 

RNsM-RM 

RM=RM/SXM 

RNsRN/SXN 

ZMI=-RM/TSM 

ZNI=-RN/TSN 

RSM=RH-HM 

RSN=RN-MN 

XN0S6=+PA0N* VON/RSN»WON/ AAN 
XMAIN=*PAOM*VOM/RSM-VOM/AAM 
THETAI=-(XN0SE*ZNI-(XMAINfZf1I))/(e+A) 
ZCGI-'XMAIN-A*THETAIVZMI 
OC  10  I=1»NFM 
QO(I)*0. 

10  C(I)=3. 

RETURN 

ENO 
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SUBROUTINE  CCEFF  (Y»  A,6.C,0) 

DIMENSION  Y<<*> 

A=Y (2) 

6=Y(1) 

C=<96.*Y<i)*56.*Y{2>-64.*Y{3)*d.*Y(4>>/(-12fc.) 
0=(-16.*Y(l) -12.*Y<2>U6.*Y<3> )/  <-12t$*> 
RETURN 
END 
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SUBROUTINE  TLCOK  (XtSLOPEtVCEPT tS«P?N) 
DIMENSION  S( 30) »P (33 ) 

THIS  IS  A  2  OIHENSIQKAL  TABLE  LOCK  UP  ROUTINE 
KITH  LINEAR  INTERPOLATION 

X  IS  THE  CURRENT  VALUE  OF  STROKE 
SLCPE  ANO  YCEPT  ARE  CALCULATED  ANO  RETURNED 
S  ANO  P  HAKE  UP  THE  TABLE 
N  JS  THE  NLHBER  OF  VALUES  IN  THE  TABLE 
CO  1  I-1«N 

IF  (X*GEoS(I)*ANO.X.LT«S(I+l))GO  TO  2 

1  CONTINUE 

2  SLOP£=CP  (I^l)-P(I)  )/(S(IU)-S(I>*.01> 
VC£PT=P(l)-SLOPE*S(I) 

RETURN 

ENO 
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PROG RAM  TAXI  <  INPUT  * OUTPUT  » TAPE 5 “INPUT *  TAPE6*0UTPUT *  7AP£3*TAPE2» 

1  T.*.PE7) 

C  THIS  IS  THE  TAX  1  HA  IN  ROUUAE  TOR  A  C-S*  AIRCRAFT  SIMULATIOA 

*****•**•'>  •#•»**»*«««•««**•«•»«•«?  •♦****  •♦**♦•♦♦•»*###•*♦•##•**♦  •#«*•#** 
♦  *#*♦•  #*#*«*•#•*«*<*#•  *•-.**>****  **#*♦♦#*»♦#*.»♦#*  *«««••»•••••**»«« 
«•«#•«*»»••*•« «**««*««V***»«*« *♦*••*#♦#*♦»***♦•**•*♦*•♦*#•*##***# *#♦***♦ 

C  ♦•♦♦THIS  OtCK  IS  FCSt  fHc  C-5A  AIRCRAFT  OAL Y**°* 

C  THIS  PROGRAM  Mil*.  PERFORM  A  SYMMETRIC 
C  DYNAMIC  TAXI  ANALYSIS  ON  A  FLEXIBLE 
C  VEHICLE  TRAVERSING  A  GIVEN  RIGID  ROANAY 
C  PROFILE. 


COMMON/F LEX*  /  SI HAIM  115)  .SIHAIN2  (15)  *  SINOSE  (15)  ♦  S  ICGdS ) 
1«SITAIL(1S?»SIPS(15) 

COPMON/F LEX2/NFM*  GH ( 15 ) *  CMEGA ( 15 ) 

C0MMGN/FLtX3/Q(15) «QG (15) .000(15) 
CGMH0N/X1/N.XM«NN,MCG.HM,MN.A,8*C*HMI 
COMMON/ X2/PA CM* VQM.AAM* AHM.OAM 
COMMON/ X 3/PA CN* VON  * A AN* AHN*OAM 
COMMON/ X A/SXM « SXN »SLP* SLA* TSM.TSA 
C0MM0N/X5/CL  *C0* AREA  * THRLST 

COMMON/ X 6/ 2 »  REACT Ml  *  RE AC IN* REACT M2  »DX» NT RUN, SPEC 0 
CCHM0N/X7/AM  *0H*CH.OH* AN.BN.CN*  OA  * AHA*  BMA.CMA.OH A 
COMMON/ X8/STR0KM  (20  « PIN CM (20) *STROKN(20 ) »PIN0A<20) *NSCM»NSCN 
COMMON/ XS/FS Mi* FSM2*  FSN* FTM1,FTN2  *FTN, XHAIN1*XMAIA2*XNQSE« 

1 VELM1* VELM2* vELN 
CR.IMOn/X  1Q/2PH,  ZPMA  ,  2PN 
CIMEXSIOA  PLANE(A) .SITEI6) 

CIMENSIOA  T ( lfl) *  T0(1A),  T00C18) 

CIHENSIOA  YP  («*>  *  YPN(A)  ,CATA  <NJ8>  .  VPA(4) 

0 IMEASIOA  ELEV(7500) 

CIMENSIOA  CGACC ( 180  0)  *  TIME(IAOO)*  PSACC(IAGO),  TIMEK15O0) 
CIPEASIOA  PRCMIBGO) *RHmFK(20) 

CIMENSIOA  SSPLOT (303).  STIMEUOO ) .OOPLOT (300) 

REAL  MCG*MH»MN,MMI 
CALL  PLOTSdOO  **0UH*3) 

PSA=0.0 

CGOUT*O.C 

MDR=0. 

x=0. 

LLL  “  1 
LL  s  0 
STOREt  =  3. 

ST0RE2  =  C. 

STCRE3  =  3. 

STORE*  "  3. 

AN  =  0 
XX  *  9 
PP  =  1.0 
M  =  1 
ITT  =  0 


C  READ  ANO  PRINT  INPUT  CATA 
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C  H= VEHICLE  HEIGHT  AT  CG  (POLNOS) 

C  A-CISTANCE  MAIN  GEAR  TO  CG  (INCHES! 

C  RsCISTANCE  NOSE  GEAR  JC  CG  (INCHES) 

C  C  =  OISTANCE  FROM  FRONT  HAIN  GEAR  TO  CG 

C  1  CORESPONOS  TO  REAR  FAIN  GEAR 

C  2  CORESPONOS  TO  FRONT  HAIN  GEAR 

C  FHUFASS  HOHENT  CF  INERTIA  (L3  IN  SEC  SQ) 

C  PLANE3  AIRPLANE  8EING  SIHULATEO  AND  GROSS  HEIGHT 
C  PSARH  *  OISTANCE  OF  PILOT  STATION  TO  CG 
C  TAILRH  *  OISTANCE  CF  TAIL  STATION  TO  CG 
C  TAKCFFs  TAKE-OFF  SPEED  (FEET/SEC) 

C  SP>.C  -INITIAL  VEL  OF  AIRPLANE 
C  THURST*  10TAL  AIRPLANE  THRUST 
C  CL=LIFT  COEFF. 

C  AREAxMING  AREA 
C  CO=CRAG  COEFF. 

C  HH* HEIGHT  OF  FAIN  GEAR  (EACH) 

C  MN=NEIGHT  OF  NOSE  GEAR 
C  SXH*  NUF8ER  OF  FAIN  GEAR  STRUTS 
C  SXN=  NUNBER  OF  NOSE  GEAR  STRUTS 
C  AHN  HYDRAULIC  PISTON  AREA  NOSE  SQ  INCHES 
C  AAN  PNEUMATIC  PISTON  AREA  NOSE  SQ  INCHES 
C  AHH  HYORAUL1C  PISTON  AREA  HAIN  SQ  INCHES 
C  AAH  PNEUMATIC  PISTON  AREA  MAIN  SQ  INCHES 
C  PAON  NOSE  STRUT  PRELCAC  PRESSURE  PSI 
C  PAOH  HAIN  STRUT  PRELCAC  PRESSURE  PSI 
C  VON  NOSE  STRUT  INITIAL  VOLUME  CU.  IN. 

C  VOH  HAIN  S « »*UT  INITIAL  VOLUME  CU.  IN. 

C  CAN  OR IF ACE  AREA  MAIN 

C  CAN  OR1FACE  AREA  NOSE 

C  SLM=MAIN  GEAR  STRUT  LENGTH  UNLOADED  INCHES 
C  DISTANCE  FROF  CL  OF  AXLE  TO  CG  LINE 
C  SLN-NOSE  GEAR  STRUT  LENGTH  UNLOADED  INCHES 
C  DISTANCE  FROF  CL  OF  AXLE  TO  CG  LINE 
C  TSM  FAIN  TIRE  SPRING  CONSTANT  PER  STRUT 

C  TSN  NOSE  TIRE  SPRING  CONSTANT  PER  STRUT 

C  OX=TIME  STEF  SIZE 

C  REAO  METERING  PIN  DESCRIPTION  STARTING  AT  ZERO  STROKE 
C  NSCN*=  OF  METERING  PIN  CHANGES  NOSE  GEAR 
C  NSCFaE  CF  METERING  PIN  CHANGES  HAIN  GEAR 
C  NFH  =  HUMBER  OF  FLEXIBLE  HOCES 
C  SIXXXX(I)  x  MODE  SHAPE  CEFLECTION  (NON  OlH.) 

C  GH(I>  =  GENERALIZED  HASS  (PGUN9S-SEC  SQ/IN) 

C  OMEGA  (I)  x  MODAL  FREQUENCIES  (RAO/SEC) 

REAO  (5*1)  PLANE 
1  FORMAT (6A1Q) 

REA0(S*5)  M«  A*  3,  C*  HHI 
5  F0RMATUF1Q.1.F12.C) 

REAO (5*10)  PSARH, TAILRH 
10  F0RHAT(2F1C.2) 

REA0(S,15)  SPEEO , THRLST , T AKOFF 
15  FORMAT (3F1Q.3) 

R£A0(5,20)  CL, AREA, CO 
20  FORMAT (3F10 • A) 

REA0(5,25)  MF,HN,SXM,SXN 
25  FORHAT(AF10.2) 

REAO(5,30)  AHN, AAN, AHH, AAH 
30  FORMAT (AF1Q ,5) 


62 


READ! 5. 35)  PACN.PAON,  VON  , VQM.OAM, OAN 
35  FORMAT (6F1Q • 5) 

READ  15*40)  SIN,SLN 
40  FORMAT  (2F1Q.  31 

R£A0(5«45)  fSH  ,TSN 
45  FORMAT (2F1Q.1) 

REAO(5*5C)DX 
READ (5*5 1) IFPLOT 
58  FCRMATCF10.3) 

R£A0(5,51)  NSCN 

51  FORMATCI5) 

F:A0(5,5?>  (STROKN(I) ,PIKON(I)« Is 1* NSCN) 

52  FORMAT  (2 (Flfl  *3) ) 

READ (5, 51) NSCN 

R£AO(5«52)  (STROKN(I) ,PIW0N(I) ,I=1,NSCM) 

REA0(5*51)  NFM 

R£AO(5,53)  (SIMAINi(I)  ,SIMAIN2 ( I } ,SINOSE  (I)  ,SICG(I)  *SITAIL(I)» 
IS IPS (  iv  »  1=1 »HFM) 

53  FORMAT (6F13 .  3> 

R£A0(5*54)  ( CHtX) t OMEGA (I) *1=1* NFM) 

54  FORMAT (FI 5* 2  »F10 .3) 

NH=NM/306 • 
fN*HN/36£  * 

MCG=M/366. 

TYPRUN  *  .01 
MTRUN=  1 

IF (SPEED  .GT.  5.)  56*57 

56  TYPRUN  =  .31 
RTRUN  *  2 

57  WRITE(6,55> 

55  FORMAT  (1H,46X, *****•***•♦•♦♦*♦**•*•♦♦•**•******♦•**•♦•**) 

WRITE (6, 60 

6G  FORMAT (46X»46K**##**#*****#**#*  INPUT  CATA  #****♦•♦#******•*) 
WRIT£(6,£1) 

61  FORMAT (///IX. 41M*******»*  GENERAL  AIRCRAFT  DATA  *•*♦**♦*•) 
liRITE  (6*  £2)  PLANE 

62  FGRMAT(//,1X,0A1Q) 

WRITE (6,65)  W,MM,NN,A,8,C»FMI 

65  FORMAT (///*5X, F10.1.5X, •««=*, F10. 2, 5X,*MN***F10. 2, 5X* 
1*A=*,FIG.3,5X,*B=*,F10.3,5X,*C=*,F10.3,5X,*MNI***F12.0) 

WRITE (6,66)  SXM,SXN,SLM,SLK*PSARF,TAILRM 
66  FOfinA T(/  /  ,5X»*SXH*#,Ftf.  1 ,5X,*3XN=* ,F9«- 1*  5X*  *SLH*#»F9.1»  5X» 

1*SLN=*,F 6 .1, 5X,*PSARM  *,F7.1,4X,*TAILRK=*,F7.1> 

WRITE (6, 7 u>  AAH* AHM, PAOM  ,  VCM  ,0AM, TSM 
78  FORMAT (// ,5X,*AA('=*,Fd»2«5X*#AHM=*,F0.2*6X*#PAOM=#*F5.2,4X* 
1*YOH**,F8.2,5X,*OAH**,F0.2,5X,*TSK  **,£10.2) 

WRITE (6, 74)  A AN, AHN*  PAON , YON, OAN ,TSN 

74  F0RMAT(//,5X,*AAN=*,Fd.2,SX,*AHN**,F6.2,6X,*FA0Ns*,F9.2,4X, 
1*YON=*,F«.2,5X,*OAN=*,F0.2,5X,*TSN  **,F10.2) 

WRITE (6, 75)  CL,CO, AREA, SF£tD, TMRLST, TAKOFF 

75  FORMAT (// ,5X »*CL=* ,F9.3,5X,*CD**,F9.3,6X,*AR£A**,Fe.2,£X, 
1*SPE£0*# ,F7. 1,4 X,* THRUST **,F7.U»  3X,*TAX0FF**,F1C .2) 

WRITE (6,6)  ISTROKM I) ,PINON(I) ,1=1, NSCN) 

76  FORMAT (//»3X»*STR0K£  NOSE  PIN  CIAHETER  •// ,2  (3X.F1C .3) ) 

WRITE (6, 77)  (STROKM(I) ,PINCM(I) , 1=1, NSCN) 

77  FORMAT!//, 3X,*STR0K£  MAI,*<  PIN  CIAHETER  *//, 2 (3X.F10 .3) ) 
WRITE  (6,60 

60  FORMAT  (1H,1X,*M00£  SIPS  StNOSE  SXCG  SIMAIN1  SIHAIN2 
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1SITAIL  CMEGA  GEM.  MASS*) 

CG  85  1=1, NFH 

85  WRITE(6»  9C  )  Cl, SIPS ( I ) , SINCSE ( I) , S ICGC I) *SIHAIN1 CX) ,SXMAINZ (I) , 
lSITAlL(I) ,QM£GA(I) *GH(I) ) 

90  FCRMAT (//,  2X,I2,7(F7.2,2X) ,F1C.1) 

c  zr.i  =initial  nain  gear  position  inches 

WRITE (6,95) 

95  FORMAT (✓// ,N6X»60H*********  INTI Al  CONDITIONS  **•*•♦•*•***) 

C  ZNI  = INITIAL  NOSE  GEAR  POSITION  INCHES 
C  ZCGI  =  IMTIAL  C.G.  POSITION  INCHES 
C  TH£TAI*INITIAL  PITCH  ANGLE  DEGREES 

C  THE  ABOVE  PARAMETERS  ARE  CALCULATED  IN  SUBROUTINE  IC 
CAc.  ICC2CC-I  .ZNI1,2MI2,ZNI,TM£TAI) 

WRITE (6, 155) 2HI1«ZHI2,ZNI,TH£TA1 , ZCGI 
105  FORMAT  (/»5X,*ZHI1=*»F7«3»5X,*ZMI2=*»F7 .3,6X{*ZNI=*»F7,3* 7X, 
1*THETAI=*,F7.3,6X,*ZCGI=*,F10.3) 

REACTN= ( 2NI* TSN  )*SXN 
REACTM1  =  CZMI1*TSM  )  *SXM/2. 

REACTH2*  (ZMI2*TSM  )  *SXM/2. 

WRITE (6, 136) XNAIM  »XMAIN2»XNOSE»REACTNl»REACTM2,RtACTN 
}  1C6  FORMAT  (/5X,*XNAIW1=*,F7 .2,3X**XMAIN2**,F7 .2 ,6X,*  XNOSE**, F7*2,5X, 

l*REACTNl=*,Flfi.O,2X.*R£ACTN2=*,F10.0,2X,*REACTN**,F10.0> 


C  READ  RUNWAY  PROFILE  DATA  (ELEV) 

C  SITE*  RU*MAY  PROFILE  ANT  DIRECTION 
C  NPTSS==  OF  NUNHAY  EL?  AT1  )N  OATA  POINTS 
WRITE  (6,111) 

111  FORMAT  CM*) 

CO  112  1=1,50 

112  £L£V(I)=0.0 
REAO (2, 1 )  SITE 
READ (2,116)  NPTSS 
NPTSS=NPTSS*50 

116  FORMAT(IS) 

L0=51 

117  LSD  *  LQ  ♦  9 

fi EAO (2 y 116)  (ELEV (I) ,I=LC,LSO) 
lie  FQRHATUCF7.3) 

MO=LO 

IF(LSO.Gt. NPTSS)  GO  TO  120 

LQ=L0+1C 

GO  TO  117 

12C  ELtVl=ELE V(51) 

CO  125  1*51,  MO 
125  ELEV(l)  =  £L£V(I)-ELEV1 
CALL  RETURNS (2) 

CISTAN=0. 

LS0*LS0-9C 
SLP= (LSO-50) *2 
SLP*ELE V (LSOJ/SLP 
CO  126  1=51  ,N10 
ELEV ( I) *£l£V (I) *SLP*OIST  AN 
126  0 1ST AN* 0ISTAN*2, 

IVAL=(A«e)/26. 

VALA  =  (A+C>/2<t. 

IVALA*VALA 
WRITE (6, 130) 

130  FORMAT (1 1*1,5 X,*  RUNWAY  PROFILE  CATA  NORMALIZED  (SLOPE  REMOVED)* 


1S0X,*F£ET  OOMN  THE  RUNHAY*) 

WRITE (6* 1)  SITE 

1.01=1 

135  LSOi  =  L01  ♦  9 

LPRIN  *  LSOI  ^  2 

WRIT£(6*14(J)  (EL£V(I)»I*L0i«LS01)  .LPRJM 
140  FORHATC3)«»10F10.5*iaX,I»> 

L01=L01«10 

TO  145 


145 

15C 


155 


ENO  RUMHAY  OAT A 


IF(LOI.GE.NIO)  GO 
GO  TO  13* 

WRITE (6*150)  LOi 
FORMAT (42H*******' 

MiT£{/»i<3IT£ 

ENORUN  =  NPT$S*2 

SlKuffiffiiSilJ^WTWiHiaoan/f.imoo. 

£NDRUN*EKORUK-90 . 

T(2)  *ZCGI 
T0(2)  *  6. 

T  (4)  *  ZM1 
TQ(4)  =  0. 

T(6)  *  ZM 
TO (61  =  c. 

T(81  *  THtTAl 
T0(8)  *  0. 

T (10 1  *  C. 

TO(10)aSPEEO 
T (12)  *  ZHI2 


(.01**19) 


T0(12)  a  0.0 
T (141*  TO  (14) 

T ( 16) *  TC (16) 
T(10)=  TC(16) 
WRITE (6» 155) 
FORMAT (iMl*7X.* 
1  TAILAC 
2£M 
cz*o. 

2PM* 0* 

ZPN*0. 

zoot=o. 

ZHO0T«0. 

ZPMA  *  0.0 
ZAOOT=0. 


*T00(14)  *  0.0 
=T00(16)  *  Q.C 
=T0D(16)  a  0.0 


XMAIN1 

SPEEO 


XKOSE 
OIST  • 


FSMi 

CGACC 


PSA 


J=1 

l&C  YP(1)  *  ZOOT 
VP(2)s£L£¥(J) 

YP(31 «ELtV( J*l) 
YP(4)*£LEV(J*2) 

YPN(l)*ZMOOT 
tPW(2)s£L£V(  IlfAL) 
YPN(3)*£L£V(ItfAL*l) 

TPH(4) *ELEV( 19AL*2) 
TPA(t)*ZAOOT 
VPA (2)  *  fcLEV(IVALA) 
yPA(3)  *  ELEVdVALAd) 

YPA(4)  *  EL6YUVALA42) 

Z*OZ 

CALL  COEFF  (YPN* AH*8H*CN  OH) 


FSN 

TIM 
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CALL  COEFF  < TP, AN,BN,CH» CM) 

CALL  COEFF (YPA.  AHA.BMA.CMA.OMA) 

Z90T*BH 

ZNOOT=BN 

ZAOOT*BPA 

IF(M.GT.l)  Rf.  245 
(n  H  ♦  i 

IBS  MRITE(6«170)  XMAINl*XNCSE»FSNltFSN»TAILAC*TO(10) «T(10) * 
1C60UT  .PSA* X 

XRITtC’.BCO)  TAILAC.CGOUT »PSA 
900  FORMAT!  >X,F1G.3) 

TIMEX* 

170  r LRHaT( 3a*£{3X«F1C.3)«2 (3X.F10.0) .5(3X.F10.3) ,2X,F6.2) 


IF(IFPLOT.EQ.l)  GO  TO  230 

IF  (A8S(CG0UT  ) .LE . ABS (STORE2) .AHC.ABS(ST0RE2) .GE .ABS( STORED ) 
1GO  TO  ISC 

175  STOREi  =  STORE2 
ST0RE2  *  CGOLT 
GO  TO  190 

ISO  IF ( X-TIME (NN) .GT . .CS)  GO  TO  185 

IFU8SCST0RE2I  »GT. .30000)  GO  TO  1S5 
GO  TO  175 
185  fch  *  NN  4  1 

CGACC(NN)  *  STORE2 

PROF (NN)  =  ZPN 

TINE  CNN)  =  X  -  10.  *  OX 

IF  (TINE  (AN)  .LT.O  • )  TIME(NN)  =  .01 

STOREI  *  STORE2 

ST0RE2  *  CGOLT 

190  IF  (ABS(PSA)  .LE.AES (ST0RE4) .  ANO.A6S  (STORE**)  .GE.A8S (STORE3) ) 

160  TO  °05 

200  STORf  »  STORE*. 

STORi  PSA 
GO  TO  215 

205  IF(X-TIM£1(LL).GT..08)  GO  TO  210 

IF  (ABS  (STORE**)  .GT..32)  GO  TO  210 
GO  TO  200 
210  LL  *  LL  ♦  l 

PSACC(LL)  =  ST0RE4 

TINE1(LL)=X  -  10*OX 

IF(TIMEKLL)  .LE.O.)  TIHEl(LL)  =  .001 

STORE3  =  ST0RE4 

ST0RE4  =  PSA 

215  IF(A8S(T(10)  -RH»  .LT.5.0)  GO  TO  223 

GO  TO  225 
22C  IX  «  II  A  1 

SHARK ( II )  =  X 
RH  *  RM  ♦  1003. 

ITT  *  II 

225  IF(ABS(X*PP) •GT..005IGO  TO  230 
SSPLCTCLLLI  =  TO <101 
STIME(LLL)  *  X 
COPLOT (LLL*  =  T (10) 

LLL  *  LLL  ♦  1 
PP  a  PP  ♦  l.fi 


230  IFCM0R.LE.80 .)  GO  TO  235 
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HOR=0. 

WRITE (6, 155) 

235  FOR  =  HOP  ♦  1. 

240  X  =  X  ♦  CX 

CALL  TAYLOR (T»T0, TDD) 

IF(T(10) .GE.ENDRUN)  GO  TC  260 
IF(X.GE.3;0.)  GO  TO  263 
TIPEXsTIPEXfOX 
QOOCG=C. 

COOPS-O. 

cootal=o. 

CO  241  1=1, NFH 

QL)OTAL=CCuT  £L+CDD(I)*SITAIL(I) 

COOCG  =  QOCCG-fQOO( I)*SICG(  I,' 

24  1  QOOPS=QCOPS*QOOCI>*SIPSCI» 

TAIL AC2 (TDD (2) *T  AILRP*TOC (8)1/ 3# 6.  ♦000TAL/386. 
PSA=  (TOO (2) -PSARP*T00(8)  )/386.  ♦Q00PS/386. 

CGOUT*TOC(2)/386.  ♦Q00CG/386. 

IFCTOCIO) .GE.TAKOFFI  GO  TO  250 
IFCZ.LT.4.)  GO  TO  245 
CZ=Z-4. 

J*J*2 

IVAL=IVAL*2 
IVALA  *  IVALA  ♦  2 
GO  TO  16C 

245  IFCTIHEX.LT. TYPRUN)  GO  TC  240 
GO  TO  165 
250  NRITE(6,255) 

255  FORMAT (5X»*  THE  VEHICLE  PAS  TAKEN  OFF*) 

260  SRITEC6.262)  ENDRUN, T (10  > 

262  FORMAT C 3X » *  ENO  OF  RUNHAV*,2F10 .3) 


TOTAL  =  X/.01 
NRITE(6,3G0)  TOTAL 

300  FORMAT (*  TOTAL  NUMBER  OF  PCINTS  CN  TAPE  IS*F12.0) 

263  HGM  *  X 

IFCIFPLOT.EQ.il  GO  TO  285 
XLCNG  z  FLOAT (HGH) 

WRITE (6, 265)  NN.LL 
265  FORMAT (2120) 

IFCNh.LE.180S.OR.LL.LE.1800)  GO  TO  266 
WRITEC6, 267) 

267  FORMAT (3X»*TPE  ARRAYS  CGACC  OR  PSACC  OR  PROF  HAVE  EXCEEOEO 
1  THEIR  OiMENSIONEO  SIZE*) 

266  CONTINUE 

C266  CALL  PLOT SC OAT A, 438) 

CALL  FACTOR  (2.0) 

CALL  PLOTCQ.0,-11. (»-3) 

CALL  PLOT (3. ,.7, *3) 

TIKE  CNN* I)  =  0.0 
TIPECNNF2)  x  l.Q 
TIHEKLL  +  1)  =0.0 

TIP£1UL*2)  *  1.0 
CGACC (NN*1)  =  -1.0 
CGACCCNN+2)  *  1.0 
PSACC (LL*1I=  -1.0 
PSACCCLL*2)  =  1.0 
CALL  SCALE  (PROF  10  •  .RN,  1 1 
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PRCF1Q  =  PR0f(NN*2)  *10. 

IF(PR0F1C.GT .10.5)  GO  TO  270 
PROF (NN*2)  =  6.0 
IF1NPTSS.GE.100C)  GO  TO  270 
00  269  1=1, NA 

269  PROF  (I) =FROF (I) *36. 

270  CALL  AXl S (0 . «0.,11HTIM£  (SEC.) »-ll,XLONG»0. 0, TIME (VN*i) * 
1TIPE<NN42),3) 

IXLONG  *  XLONG 

CALL  PLOT  (XIONG, 1.1  m 

CALL  PLOT  (0. .1.1*2) 

r ALL  SYMBOL  ( . 1# .4* . 105, OHOISTANCE , 0 . * 6) 

CALL  SYMBOL  (.1, .2, . 105, 5HSPEE0, C. ,5) 

CO  275  1=1, IXLONG 

CALL  NUMBER ( ST IME  (I)  ,  .2*  .  105, SSPLOT ( I)  ,0  .0, 4HF4. 0) 

27 5  CALL  NUMBER (STIPE (1) , *4, .105,OOPLOT (I) , 0  «Q,4HF5. 0) 

CALL  PLOT  (XLONG, 1.5,3) 

CALL  PLOT  (0  .,1.5,2) 

CO  280  1=1, ITT 

280  CALL  SYM£OL(fiHARK(I) ,-.l, .245, 70,0. 0,-1) 

CALL  PLOT  (XLONG, 1.9, 3) 

CALL  PLOT  (0  ., 1*9,2) 

CALL  SYMBOL  (A  .,9  .,  ,1*1,  PLANE,  0.«4C) 

CALL  SYMEQL(4.,8.,.1«»,SITE,0.,4Q) 

XL0NG2=XL0NG/2. 

CALL  PLOT(0.,1.5,~3) 

CALL  AXIS  (0,  ,~1. *3HCGA*3«2. 0,90  .  ,CGACC (NN*1) ,CGACC(NN+2)  ,*1) 
CALL  NUM6ER(0.4«-.4,.105,-.4,00.0,4HF4.1) 

CALL  NUMEER(0.5,*.4,.105,4.4,00.0,4HF4.1> 

CALL  PLOT(0.,*1.0,>3) 

CALL  LINE (TIME,CGACC,NN«1,0«74) 

CALL  PLOT  (XLONG, 3. 1,3) 

CALL  PLOT  (0 .,3.1,2) 

CALL  PLOT(0.,3.5,~3) 

CALL  AXIS(0.,<*1.,  3HPSA»3 ,2.0 ,90  .*  PS  ACC  (LL*1)  «PSACC(LL+2)  ,-l) 
CALL  PLOT(0. ,0.,3) 

CALL  PLOT (XLCNG*Q.«2) 

CALL  PLOT  (ALONG, .4,3) 

CALL  PLOT  (-•!•• 4, 2) 

CALL  NUHB£R(0.5«  .4,  .105, *.4, 00.0  ,4HF*».l) 

CALL  NUH8ER( Q.4*-»5» . 105 ,-.4, 00 • 0,4HF4 .1) 

CALL  PLOT(Q.,-1.0,-3) 

CALL  LINE (TIP£1,PSACC,LL,1,0,  0) 

CALL  PLOT(0.,-3.0,-3> 

CALL  AXIS(-t.5,0.,l5HELEYATION  (IN.), 15, IQ. 0,90.0, PROF CNK+1) , 
1PR0F(NN*2),-1> 

CALL  PLOT (Q««Q»*~3) 

XPR0F*.25+ABS (PROF (NN*1) )/PROF (NN+2) 

CALL  SYMBOL  (2., XPROF, ,14,15MN0SE  GEAR  TRACK*0.*i5> 

CALL  LINE (TIME, PROF, NN, 1,0,  0) 

>ST0P*ALCNG45. 

CALL  PLOT  (XSTOP,Q.,-3) 

NRITE (6 ,2 90) 

CALL  PLOT (0.,0.,40) 

C  CALL  PLOTE 
285  STOP 
290  FORMAT(*S*) 

£NO 
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SUSROUTIhc  T AYLOR (T»TQ, TOO) 

COKMCN/FLEX1/SIMAIN1 (15) ,SIMAIN2 (15) ,SIN0S£(15) ,SICG(15) 
1,SITAIL(15), SIPS (IS) 

C0MMCN/FLEX2/NFM ♦ GH ( 15) , OMEGA ( 15 ) 
COMMON/FL£X3/Q(15),QO(15),GDO(15) 
C0MM0N/Xl/M,)iM,NN,HCG,NH,NN,A,8,C,NNI 
COMMON/ X 2/PA CM, VOM  ,AAM,AHM,OAM 
C0MM0N/X3/PACN, VON, AAN, AHN,OAN 

C0MM0N/X6/Z, REACTH1 ,R£ACTN«REACTM2,QX,NTRUN*SPE£C 
COMMON/ X5/CL ,C0, AREA, THRUST 
COMMON/ XG/SXM,SXN,SLN,5LN,TSM,TSN 
CCMMCN/X7/AM ,8M,CM  «QM ,AN ,8N ,CN, ON « AHA, 8MA«CMA,0M A 
Cji1M0N/xe/ST90KM(2C> ,PIN0M(2G) ,STROKN(20) ,PINON(20) ,NSCM,NSCN 
C0MMQN/XS/FSM1,FSH2*FSN,FTM1«FTM2*FTN«XMAIN1,XMA1N2,XNQSE, 
lVELMl, VELM2, VELN 
COMMON/ X15/Z PM, ZPMA,2PN 
CIMEKSIOF  T(18>,  TO ( 16) ,  TC0(18) 

REAL  MCG , MM, MN, MMI 


C  THIS  THE  TAYLOR  SUBROUTINE  FOR  A  C-5A  AIRCRAFT  SIMULATION 


1  ZsZdQ(lC)*0X+TQQ(13)*0X**2/2. 

2PN=AN+8N*Z+CN*Z**2>0N*Z**3 
ZPM=AM48M*Z^CM*Z**2+0M*Z**3 

ZPNA=AMA*8MA*Z+CMA*Z#*2*0MA*7**3 
QTN=0 • 

QTM1=  0.0 
CTM2=Q. 0 
CTON=0. 

CTOM1=0.0 
CTOH2=0.0 
00  10  1=1, NFM 
QTN=QTN«0(I)*SIN0SE(I) 

CTH1=  QTM1  ♦CCH  *S  IMA  INI  (I) 

CTM2=  QT M2  ♦ C ( I ) *S IHAIN2 (I) 

0T0M1*  QT0M1  ♦  00(1)*  SIMAINl(I) 

CTCM2=  QT0M2  ♦  00(11*  SIMAINZd) 

10  QT0N=QTCN*QC(I) •SINOSE (I) 

XNCSE  =  (T(21  -  3  *  T(8)  -  T(6)l  ♦QTN 
XMAINis  (T (2 )  ♦  A*T(5)-T(4))  ♦  QTM1 
XMAIN2=  (T(2>  -  C*T(8) -T (12) )  ♦0TM2 
VELN  *  TC(2)  -  B  *  T0(8)  -  TQ(6)  ♦  QTON 
VfcLMl  =  T0(2I*  A*T0(8)-T0(«»)  ♦QT0M1 

VELM2  =  T*3(2)-  C*T0(6)-T0(12)  ♦CT0H2 

CALCLLATE  SPRING  ANO  OAMPING  COEFFICIENTS 

IF (XNOSE *GT . C .)  XNOSE  =  -.001 
IF(XKAIN1*GT  .0.)  XMAIN1  =  -.001 
IF(XMAIN2.GT.O.)  XMAIN2  =  -.031 
IF ( VELN. EQ.O  .)  VELN=,001 
IF(VELMl.tQ. .00)  VELM1*. C01 
IF(VEL^2.EQ.C.)  VELM2=C. 001 
VELN  = VELN  -  TO( 14) 

VtLMi“VELMl-  TO( 18) 

V€LM2*VELN2-  T0U6) 

XNLK=A8S (XNOSE) 

XHLK1  *  A3S ( XHAIK1) 


70 


XMLK2  *  AQS(XHAIX2) 

CALL  TLOCMXNLK.SLOPeN,YCEPN,STRCKN»P!NON,hSCN) 

CALL  TLOCKCXPLXlt  SL0P£H1,YC£PH1*£TP0KH«;  THOHtNSCr!) 

CALL  TLQCK(X£LK2 « SL0FEM2  * YCEPN2 » STRQKHtPINuM.NSCN) 

AON  =  SLCPEN* XNLK  ♦YCEPN 
AOM1  =  SL0PEP1*XKLK1  ♦YCEPK1 
ACH2  *  SL0PEF2*XHLK2  ♦YCEPH2 

CONa(,0Q0Q8* (AHN**3.)  )  / <2.* C .9  * >ON#**2) 

COMl=  (.03008*  (AFH**3.) ) / (2  »*( .9*A0H1 )  **2) 

C0M2=  (.QC00fl*(AHM**3.))/(2.*(.9*AGM2)**2) 

SSN=(PAON*VfOM/(  ( C  VON/AAN) -XNLSO  ) 

SSK1=  (PAOH*VOM*/ ( ( ( VOH/ AAH) -XHL*1> ) 

SSM2s  (PAUM*  yOMI / (  (  { VOK/ AAH) -XNLK2! > 

FTN  =  SXN*  TSN  •  (T(6)  -  ZPN) 

FTM1  *  (SXH/2.)*TSH*(T(4)-ZPM) 

FTM2  -  (SXN/2.)*  TSM* (T(12) -ZPNA) 

IF (FTN.GT .0. ) FTN=Q  • 

IFCFTH1.GT.0  .1  FTH1=0.0 
IFCFTH2.GT.0 .)  FTH2=0.0 

IF(T(14) .LT.C.000)  SSN  =  (2937242 . )/ i i 13 .35) -XNLK) 

IF(T(16) .LT. 0.900)  SSM2  =  (5124326. 59) /( (37.893) -XMLK2) 

IF (T (18)  .LT.C.000)  SSH1  =  (512432G .96) / < (37 .89) - XHL Kl) 

FSN  *  SXN* (-SSN*CON*VELN*ABS  (VELM ) 

FSKi  =  ( S XM  /2.C)  *  (-SSH1  ♦  COMl*Y£LNi*ABS WELN1) ) 

FSH2  a  (EXH  /2.0)  *  (-SSH2  ♦  COH2*VELH2* A8S (v£LH2>) 

VLIFT  =.CC1189*CL*AR£A*(TD(10)*TD(1C)) 

CR AGA* VL IFT*CO/CL 

CRAGT  =A6S(.025*FTM1  ♦.025*FTH2  ♦.G25*FTN) 

IF(NTRUN  .EQ.  1)  GO  TO  125 
IF(T0(1C)  .LT.  SPEED)  GO  TO  125 
THRUST  a  ORAGA  ♦  ORAGT 
C  SECONDARY  PISTON  CALCULATIONS 

125  CON  =  .5 
COH1  *  .5 
COM2  =  .5 

IF('f 0(14). GT. 0.0)  CON  a  S.C 
IF (TO (16)  ,GT .0*0)  COM2  *  5.C 
IF(T0(18)  .ST.O.Q)  CONI  a  5.0 

F2N  a  2733218. / (18.43-AES(T (14 ) ) )  -20 ,*TO(14) -CCN*T0(14)*ABS( 
1T0(14) ) 

F2H2  a  4E43757./ (26.21-AES(T(16) ) )  -20 .*TO( 16) -CCH2*TQ (16) *A8S( 
1T0(16) ) 

*"2H1  a  4E43757. /  (26 .21-ABS (T  (18 )  ) )  -  2G.*TOC18)-COM1*TO(18)*A8S( 
lTO(lfl) ) 

FSTN  a  fsn*  .803  ♦  F2N 

FST1  =  ( FSH1/2.) *• 708  ♦  F2MI 

FST2  a  (FSH2/2.)  *.738  ♦  F2H2 

T00(2)a(-FSN-FSN1-FSN2-NCG*386.*VLIFTI/MCG 

T00C4)  a  (FSF1-FTH1-HN*772.)  /  ( NM*2 . ) 

TOO (6)  =  (FSN-FT  N-HN*386 • )/HN 

TOO (8)  a-CFSPl*A  -FS«2*C  -FSN*B  -0RA6T* (SLM+XHAINl) I/NHI 
TOO ( 10)  =  (TFRUST-ORAGA-CRAGT  ) /( (HCG*NH*HNJ  *12.) 

TOO ( 12) =  (FSF2-F  TN2-KH*7  72 . ) / ( HM*2 . ) 

TOO (14)  a  FSTN/. 259 
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TOO (16)  =  FST2/.259 
T00(18)  =  FST1/.259 
CO  20  I*1»NFM 

CDO(I)  =  -(SIfcAINl(I)MFSH-REAr.TKl)  ♦SlHAIN2(I)MFSH2-REACTM2) 

14  (FSX-REACTfO^SINOSE  (I)  ♦ «l*OM£GA (I)*GN ( I ) #Q0( I) 4CMEGA ( i)**Z*Z* 
2GM(I)#Q(I))/£H(I) 

OX  =  .001 

IF(FSTN.tT.8000..0R.T(1.4i  .LT.Q  .0)  OX  £  .00005 
IF(FSTi.LT.8000..OR,T(18).LT.0 ,0)  OX  £  .00005 
IF1FST2 .LT.3000.«OR.T(lb).LT.O.O)  OX  £  .00005 
IF(A8S1TC(14)).GT..001)  CX  £  .9000? 

IF(A9S(TC(16l).GT..a01)  OX  £  *00005 
IFIA03(TC<16>;  .GT..001»  CX  =  .  00005 
CC  30  1=2,18.2 

T  ( I)  =  T  < II  ♦  TOU)*OX  *  (TQQ(I)#DX**2S/2» 

TO CI>  £  TU(I)  ♦  TOO(I)*OX 

IF  (T  (141  .G£.*.OG31.ANO.,TC(14),GE.:.OQQ>  T  ( 14)  =TC  ( 14)  =  T00  (14)  =0 .0 
IF CT (16) 00 Ql. AND. TC 116). &E.;. 000)  T (16) *TO C16)=T00 (16) =0.3 
IF (T (18)  .F,E.~.00  01.ANr/.TC(18)  .GE.C.COO)  T(18)=T0(18)=T00(18)  =0.0 
CO  43  1=1, NFK 

C(t)=Q(I)4QO(i)*OX4(f/DO(I>*OX**2>/2. 

CO(I)=QO (I) +QOO( l)*OX 

RETURN 

END 
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SUBROUTINE  IC (ZCGI , ZMI1 , ZMI2, ZNI  .THETAI > 

CCMMCN/X1/H, NH,HN,MCG»MM  »MN» A»B»C ,MMI 
C0MMCN/X2/PACM, VQM,AAM« AHM.OAH 
C0MM0N/X3/PACN, VON, A AN* A^NyOAN 
CGMM0N/X4/SXM  VSXN, SLKySLNyTSH,  TSN 
C0MM0N/X6/Z, REACTM1 , REACTN, REACTM2 » OX, NTRUN 
C0MM0N/X9/FSHiyFSH2  *  FSN»FTHi,FTM2yFTNy  XtlAINl,  XKAIN2  yXNOSEy 
1 VELH1, VELM2, VELN 


C  THIS  IS  THE  IC  SUBROUTINE  FCR  A  C-5A  AIRCRAFT  SIMULATION 


C  THIS  PROGRAM  HILL  FIND  THE  INITIAL  CONDITIONS  FOR  TAXI 
C  FOR  SXM  MAIN  ANO  SXN  NOSE  GEAR 
C  ZMI*  MAIN  GEAR  TIRE  OEFLECTION  35  PERCENT 
C  ZNI-  NOSE  GEAR  TIRE  OEFLECTION  35  PERCENT  OEFLECTION 
C  ZCGI*  CG  OEFLECTION 
C  TH£TAI=  PITCH  ANGLE 

C  XMAIN=  MAIN  LANOING  GEAR  STATIC  STROKE 
C  XNOSE*  NOSE  GEAR  STATIC  STROKE 


1  Kfll  =  600C0. 

10  RM1=RM1  450. 

CM2  *  (RM*(A+8)  -  N*8)/(C-B> 

RN  *  H-RM  -RM2 
RM1T  *  RM1/2. 

RM2T  =  RN  2/2 * 

ZNI*  -RN/TSN 
ZMI1  *  -RM1T/TSM 
ZHI2  *  -RM2T/TSM 
RM1T  =  RM1T  -  MM 
RM2T  *  RM2T-MH 
RN=RN-KN 

XNCSE*PACN*VCN/RN  -  VON/AAN 
XMAI(a=PAOM*VOH/RHlT-VOM/AAM 

IF (RM1T.GT. 252000.)  XMAIN1  *  5223296. /RM1T  -  37.89 
IF (RM2T.CT. 252000.)  XMAIN2  *  5223296. /RM2T  -  37.89 
THETAI* (  CXMA IN1+ZMI1) -(XNOSE+ZNI ) ) / (B+A) 

XMAIN2T  =  XMAIN1-THETAI* (A+C)  -ZMI2  4ZMI1 
ZCGI  *  XMAIN1-A*THETAI4  ZMIl 
IFIRM1.GT .450000.)  GO  TO  20 
IF (A6S(XMAIN2-XMAIN2T) .LT..006)  GO  TO  3C 
GO  TO  10 

20  NRITE(6,6Q) 

40  FORMATI1X ,*N0  INITIAL  CONDITIONS  HAVE  BEEN  FOUND*) 

33  IFCZNI.GT .0.0)  GO  TO  10 

RETURN 
ENO 
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SUBRCUTHc  CCEFF  CY,  A«6tC«0) 

CIMENSIOf  V <-41 

*********************************************  **************  **  ******** 

C  THIS  IS  THE  CQEFF  SUBROUTINE  FOR  A  C-5A  AIRCRAFT  SIMULATION 

****************************** ************** ******* ****** ************ 

1  A=Y<2) 

e*Y<n 

C=(?E.*Y(l)  +  56.,‘  r(2)-64.*Y<3)  ♦8.*Y(4))/(-12d.f 
C=<-16.*Y<1> «12*#Y  C2 ) *16 •*Yt3l»4»*Y{4) )/<-12o.) 

RETURN 

END 
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SUBROUTINE  TLOOK  (X, SLOPE, YCEPT,S,P,N> 

CIMENSION  S(30),P(3Q) 

C  1HIS  IS  THE  TLOOK  SUBROUTINE  FOR  A  C-5A  AIRCRAFT  SIMULATION 


THIS  IS  A  2  GIHENSIONAL  TABLE  LOCK  UP  ROUTINE 
fcITH  LINEAk  INTERPOLATION 

X  IS  THE  CURRENT  VALUE  OF  STROKE 
SLOPE  ANO  VCEPT  ARE  CALCULATED  AND  RtTURNEO 
S  JNO  P  HAKE  UP  THE  TABLE 
N  IS  THE  NUMBER  CF  VALUES  IN  THE  TABLE 
CO  1  1=1, N 

IF<X.ftE.S<I) .AND.X.LT.S(I*l))GO  TO  2 

1  CONTINUE 

2  SLCPE=(Pim>-PCin/<sa*i>-sm*.ai> 
YCEPT=PU)-SLOPE*SCI) 

RETURN 

ENO 
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SUBROUTINE  TAYLOR (T ♦  TO.TCO) 


C  THIS  IS  TH£  TAtLOR  SUBROUTINE  FOR  *  F-4  AIRCRAFT  SIMULATION 

COMMCN/FLcXl/ SIMA  IMIS)  ,  S INOStt 1 5  )  .SICG(15)  ,SITAILI15>  ,SIPS(l5> 
CCMMGN/FL£X2/NFM,GM(15),CH£GA(15) 

COMMON/ FLEX3/QC 15) «Q0(15) «QUOC 15) 

CCMMGN/Xl/W.NM, NN ,MCG , MM ,MN . A , 9 , MM  I 

COMMON/ X2/PACH,VOH,AAH,AKH,OAM 

COMMON/ X3/PACN, VON, AAN,AHN»OAN 

COMMON/ XR/SAM,SXN, SLH.SLN, TSH  ,TSN 

I'lMM'-N/XS/CL  ,CO,A«£A,THRLST 

COMMON/ X6/Z, R6ACTH, REACT N»CX«NTRLN, SPEcO 

COMMON/ X  7/AM ,8M«CM»0M» AN. BN.CN, ON 

COMMON/ X 6/ST ROKM ( 20 ) .PIN  CM (20  J , S t  RGKNC  2C ),PINON<20) ,NSCM,NSCK 

common/xs/fsn,fsk,ftm*ftk,xmain,xkos- , V  ..1,VELN 

COMMON/Xi;/ZPM,2PN 

CIMEKSION  T ( 12) ,  TO ( 12) *  TCOCL2) 

R£Al  MCG ,MM,MN,MMI 
1  Z  =  Z»IOUC)*CX*TOOtl3)*OX**2/2. 

2PM=AH»BM*Z*CM*Z#*2*CM*Z4*3 

ZPN*AN*BK*Z*CN4Z442»ON4Z443 

QTN=0. 

CTM=3. 

CTON=iJ. 

ctom=o. 

OC  13  1=1, NFM 

QTN=QIN4QCI)*SIK0SECI) 

CTM  =  £)MC(I)»SXMAIMI> 

QT0M=CTCM*QUI)*SIMAIN(D 
10  QTOK=QTCN*QCCI)*SINOS£U) 

XNCSC  =  (1(2)  *  8  *  T(d)  -  T (6) )  *QTN 
XM.AIN  =  (T  (2 1  ♦  A  •  f((t)  -  T(4,l  tOTM 
VELM  =  T  C  (2)  ♦  A  *  TOCO)  -  TO(4)  *QTOM 
VELN  =  T C (2)  -  9  •  TO (A)  -  TOCS)  ♦QTON 
IF(XMAIN.Gc. 0.)  XMAIN=-.l 
IF(XN05£.6S. C.)  XNOS£=-.l 
IF(V£LM.£Q.O .)  V£LM=-,1 
IF(V-LK.£Q,C.)  VSlN=-.1 
XMLK=A«S(XMAIN) 

XKLKsABS  (XNOSt) 

C  NOSE  AND  MAIN  DAMPING  CO£FF 

CALL  TLOC< ( XMLK, SLOP£M, YCcPM , STRCRM, PIKOM, NSCH) 

CALL  TLOCX(XKLK,SLOP£N,YC£PN,STRCKN,PINON,NSCN) 

ACM  =  SLCPt«4A9S(XMAlN»  ♦  YCtPM 
A'JN  ?  SLCPtN**  ' XNOSt)  ♦  YE'PN 

IF  (A53(XMAIN)  ,  .  .  13 . 96 .  A.\0  .  VtL  M  *  L  T  .  0  .0  )  AOM  *  ACM  ♦  .04906 
IF  (VSLN.GT.0.0)  AON  =.6C52-AON 
IFCVtLM.GT.C.3)  AOM  =  .655  -  AOM  ♦  .4536 
CUN” C . 0  0  C  j»* CAMN*»3.))/C2.*C.9*AGN)**2) 

COM* (. Q C U6*  ( AHM443. >  )/(£.•(. 9* AU)**2 1 

IF  (V£LM.CT.O.O.ANO.A0S(XMAIN)  .GT.13.66)  COM.  =  COM  ♦  16.72 
C  NOSt  ANO  MAIN  STRUT  PNELMATJC  FORCES 

IF(»9S(XMAIN)  .GT.13.90)  SSM  =  ISSsi.  73/(2. 292S-A8SUMAINU3. 06) ) 
IF  (AgS(XfAIN) ,LT .13.66)  SSM  =  |PACM*VO«) /( C ( WOMV AAH) -A6S (XMA IN) ) ) 
IF(A9SCXMAIN) .GE.13.66.AK0.ABSCXMAIN) .LS.13.93)  SSM  =  (1865. 
i ♦  17l0aS.*(A6S(XMAIN)-13.46)) 

SSN=(PAOK4VQK)/( ( (VON/AAM  -XNLK) ) 
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FTH  «  SXF  •  7SM  *  <T(4)  -  ZPH) 

FTN  =SXN*  TSF  ♦  <T<6>  -  ZPN) 

IF(FTH.GT.Q.)FTH=Q. 

IFCFTN.GJ.0.)FTN»3. 

FSN=SXN*  (  -SSMCCN*  VELN*ABS(  VELN)  ) 

FSH=SXNM-SSN  ♦C0H»VELH*A8S<VSl.H>  ) 

IF  (AeSUFAIH)  .GE  .15.881  FSN  *  FS*  *100  00003  .♦  (XP*IN  ♦  15.58) 
IF (FSH.GT .0.0)  FSN  =  0.0 
IF(FSN.GT.O.C)  FSN  =  0.0 
VLIFT  =.OCllg9*CL*AREA*<TOtlQ)*TOI10)> 

ORAGA= VL IFT*CO/CL 

CRAGT  =A£S{.C25»FTH*.025  ♦FTN  ) 

TCCNTRUN.EQ. 1)  GO  TO  20 

IF CTOCIO )  .LT •  SPEED)  GG  TG  20 

THRU$T=DRAGA*ORAGT 

23  TOG (2)  =  <-FSN-FSH-NCG*3e6.+VLIFT)/HCG 
T  DO  f  <*)  =  CFSN-FTH-386.*SXH*Hf')/<HM«SXM) 
TOOC6)=CFSN-FTfc-NN*356.*SXN)/ t«N*SXN> 

TOO (8)  =-<FSN#A  -FSN*3  -CRAGT  • {SLH*XHAIN))/HMI 
TOO  1 19)  =  (TfcPUST-ORAGA-CRAGT) / C  INCG*12« ) ) 

CO  30  I*?.»NFH 

30  QDO<I)  =  -<SIFAIMI)MFSK  2EACTH) ♦SINOSE C I ) * (FSN-REACTN) 

1+ .10*OKECA(I  )*QO  (II*Grt(II ♦CNEGA(I)**2*GH  (I)*Q(I))/Gh<I) 

CX= . 301 

IF(A0S(Xf'AIN).GE.13.85iAAO.ABSCXFAIN)  .LE.13.91)  CX  s  .  CC  01 
CO  40  I  =  2.10,2 

TCI)  *  TCI)  ♦  TO<I)*OA  ♦  CT00tI)*0X**2)/2. 

43  T II « I )  =  TOJI)  ♦  TQD(I)*OX 

CO  50  1=1 , FFN 

Cei)=QCl)^Q0CI)*?X*C000(I)*0X**2)/2. 

50  CO  Cl >  =QO (I)+QOO(I)*DX 

RETURN 
END 
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SU9R0UTIKE  ICC  Zv't'ltZMltZNI'THElAI) 

COMHCN/F  LcXl/SIH A  IN  C 15)  *  S INOSc  C 15 ) «  S iGG ( 15) , SIT AIL C 15) « SIPS  CIS) 
C0KM0N/FUX*/NFfc,GM(15>  ♦CHEGAC15J 
CCMMCN/FL'-X3/QC15),QOC15),ODOC15) 

COMMON/ Xl/N,  KM,  NNfKCi.  ,MM,MN,A,  8,  KM  I 
C0MM0N/X2/PACM, VOM • A AM, AHM*OAH 
C0MHCN/X3/PAGN, VON*  *A-', lt-N*OAN 
COMM.ON/XN/SXM,SXN»  SL.i:  SEN  ,*TSN  ,TSN 
COMMON/ X9/FSM«FSN* FI  FI  J  .  >. MAIN,  XNOSc,  VEIN,  VELN 

C  THIS  IS  THE  IC  SUBROUTINE  Kc  A  F-4  AIRCRAFT  SIMULATION 


C  THIS  PROGRAM  MILL  FIKO  THE  IKITIAL  CCNOI'IONS  FOR  TAXI 
C  FOR  SXM  MAIN  ANC  SXN  NOSE  GEAR 
C  ZMI  =  MAIN  GEAR  TIRE  DEFLECTION  35  PERCENT 
C  ZNI=  NOSE  GEAR  TIRE  OEFLECT 1CK  35  PERCtNT  DEFLECTICK 
C  ZCG 1=  CG  OEFLEC  TION 
C  THt T A  1=  PITCH  ANGLE 

C  XMAIN=  MAIN  LANDING  GEAR  STATIC  STROKE 
C  XNOSE=  NOSE  GEAR  STATIC  STROKE 


1  RM-H/C1.4A/8} 

RN*H-RM 

RM-RM/SXM 

RN=P.N/SXN 

2MI=-RM/TSH 

ZNI*-RN/TSN 

RSH=RM-MK 

SSNsKN-M* 

XNOS£=*P AON*  MON/RSN- VON/ AAN 
XMAIN*aPACM*KOH/RSM-VOH/AAM 
lFCABStXMAlNI.LT.13.ed)  GO  TO  10 
XHAIN  x  a. 

20  XM AIK  =  AMAIN  -  .001 

RSMT  x  CPACM* VOH) /C  C  CVOH/AAH) -*ABS(XHAIK!)) 
IFCA6SCXMAIN). GE.13.ee)  RSMT  x  ♦  C195K6.)/ 

1  ( ( 1 2.298 )  -ABSCXMAIN  ♦  13.65))) 

IFCA8SCRSH  -  RSMT )  ,L T .50 . 0 )  GO  TC  1C 
GO  TO  20 

10  THSTAI=-CXNOSE*ZNI-<XHAIMZMI)  I/C9+A) 

ZCGI=XMAIN-A*THETA  !♦:«! 

DC  30  1=1, KFH 
QO(I)=0  . 

30  C(I)=I). 

RETURN 

END 
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SUOROUTm  TAYLCR(T*TO*TCOI 

CGKMCN/F  LlXI/SIK AI M  15  i » S I NOSE  (15)  • SiCG ( 15)  ,SITAIL(15>  .SIPS (15) 
CCKKCN/FLcX2/NFK»GM(15).CHE&A(l5) 

C0KKCN/FLwX3/Q(15)»QC(15) .000(15) 
COMMON/X1/W,NM,MN,KCG,MM,MN,A,0,KMI 
COKHCN/X2/PACM. VON . AAM.  AKH.  OAM 
CCKK0N/X3/PACN. VON. AAN. AFN.OAN 
COMMON/ X4/SXK.SXN.SLM.SLN.TSM  .  TSN 
COMMON/ X5/CL  ,C0«  AREA  , TMF.LST 
COKMCN/Xfc/Z. REACTH. REACTN.CX.NTRCN 
COMMON/ X7/AM *8M«CM«0M* AN «8N,CN«0K 

COMMON/XC/STFOKK (20) ,PINGH(20) . STROKN ( 20 ) » PINON (20) .NSCM.NSCN 
CC‘*MON/  XS/FSF  »FSN»  FTK.FTN.  XMAIN.  XNOSE. VELM.  VELN 


COMKCN/XIO/ZFN.ZPN 
COMMON/ Xl 1/S SN.O TDM.  CTK 
DIMENSION  T(12>.  T0(12>,  T00U2) 
REAL  MCG  .HH.MN.KMI 
C  F-lll  TAYLC*  SUBROUTINE 
1  Z=Z*T0(lC)*0X*T00(13l*QX**2/2. 

120  ZPK=AH»8K»Z*CN*Z**2*0M*Z*,*3 

2PN=AN*8N#Z+CN*Z*#2*0N*Z##3 
QTN=0. 

CTM=0. 

CTON=0  * 

QTOM*0. 

00  130  I=1«6FM 
QTN=QTN«Q(I)*SINOSE(Z) 
CTK=QTM4C(I)  •SIKAIMI) 

QTOM*QT GK+QO ( II *S IMA1N  < II 
130  QTON*QTCN+QC (I ) *5INCS£ ( II 

XNOSE  =  (T(21  -  9  *  T (61  -  T(6I> 
XMAIN  *((T(2)  ♦  A  *  T (61  -  TC4II 
VELN  =(TC(21  ♦  A  *  TO (81  -  T0C4) 
VclXi  =  T C (21  -  B  ♦  TO(8l  -  T0(6) 
IF4  ..  AIN.GE.  £.1  XMAIN=-.l 
IF (XNOSE  «GE«  0  *)  XNOSE*-. 1 
IF(v£LM.EQ.O.i  VELM*-.l 
IF ( VELN • EQ.O  •  I  VELN=-.l 
XMLK=ABS (XMA1NI 
XNLK=A8S(XN0SEI 


♦qtn 

♦OTM  I*. 66 

♦  QT OM) *  *66 

♦  QTON 


G  NOSE  ANO  MAIN  OAHPING  CCEFF 

CALL  TLCCMXKLK.SLCPSH.YCEPN.STRCKM.FINOM.NSC  f> 
CALL  TLOCMXNLK.  SLQPEN. YCEPN. STRC  KN* PI  NON* NS ONI 
AOM  =  (SLCPEK*XMLKFYCEPM| 

AON  =(SLCPEN*XNLK»  YCEPNI 
C0N=(.CCC:8* (AHN**3.))/(2.M.9*ACN>**2) 
CGM=(.00C:8* (AHM**3.)I/(2.*(.9*ACM)**2) 


C  NOSE  ANO  MAIN  STRUT  PNEUMATIC  FORCES 

SSK=(PACF*VOM/(  ((VCM./AAM1-(XNLK-14.475)  >1 


SSN=(PACN*VOM/(  ((VOX/AAM-XNLKII 

IF (XMLK.GT .0 • .AJO.XKLK.LT ^lv.  275) SSH=4790 .*201 .7*XHLK 
IF  CXS'LK.Gi.l4.275.AN0.XMUC.Lc.l4.4  75ISSH*7669.  +  lc944u. 
1* UMLK-14.275) 

fTM  =  SXM  •  TSM  *  (T (41  -  ZPMI 
FIN  =SXN*  TSN  •  (T (6 1  -  ZPNI 
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IFCFTM.GT.C.)FTH=0. 

IFCFTN.GT.0.)FTN=9. 

FSN=SXNM-SSMCON*VELN*AeSCVtLN>  ) 

FSH=SXH*  C-SSI*  ♦CON*VELH*ABSCVcLN)  )*.66 
WLIFT  =  .C01ia9*CL*AREAMTOCl'j)*TCClC)) 

CfiAGA*VLIFT*CO/CL 

CRAGT  =AESC.025*FTN».C25*FTN) 

IFCHTRUN.EQ. 1)  GO  TO  125 
TH«UST=ORAGA*ORAGT 

125  TC0C2)  =  C-FSN-FSH-HCG*30<j.>VLIFT)/HCG 
T00C4)=CFSM-FTM-336.*SXK*HH)/CHH*SXH) 
T00C6)=CFSN-FTN-HN*386.*SXN)/CHN«SXN> 

TOO (3)  =  -  CFSF*A  -FSN*8  -CRAGT  *CSLH*XHAIN) )/HMI 
100**0)  =  C  D'R’JST -OR  AGA- CRAGT )  /  C  CFCG*12. ) ) 

CO  200  I-l.NFN 

2  C  0  QOO C I)=- CSIFAIN  CD* CFSH-REACTM) ♦SINCScCI)* CFSN-REACTN) 

1 ♦ »1C#0MEGA{ I) *QC C I)*GMC I ) ♦OMEGA C1)**2*GHCI)#QCI) )/GHCI) 

cx=.aoi 

IFCXFLK.GE.14.275.AN0.XHLK.LT.14.475)  OX=.OOC1 

co  loai  i  -  2,10,2 

TCI)  =  TCI)  ♦  TOC  I) *CX  ♦  C  T  00  Cl) *DX#*2) /2. 

1001  TOCI)  =  TOCI)  ♦  TOOCI)*OX 
CO  1002  1=1. XFH 

C  CI)=QCI)  +QO  CD*OX+CGOOC  l)*0x**2)/2. 

1002  COCI)sQOC I)  *GOOC I) *OX 
RETURN 

ENC 
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SUBROUTINE  IC(  ZCGI»ZMI,ZM,THETAIJ 

COMHON/F IEX1/SIMAIM15) , SI NOSE (15) ,S2CG(15) ,SITAIL<15> ,SIPS<15) 
C0MM0N/FLEX2/NFH , GM ( 15)  , OMEGA <155 
COMMON/ FI EX3/Q(15) ,Q0(15),Q 0D<15) 

COMMON/ X  1/M,  NM,NN,MCG,MM,NN,A, B,MMI 

COMMON/ X2/ PA  CM «VOH,AAM, ARM ,OAH 

COMMON/X  3/PACN ,VON , AAN, AHN, OAN 

COMMON/ X4/SXM • SXN , SLM, SLN ,TSK  ,TSN 

COMMON/ XS/FS7,FSX*FTM,FTN,XMAIN,X NOSE, VELM, VELN 


C  F-lll  IC  SUBROUTINE 

C  THIS  PROGRAM  MILL  FINO  THE  INITIAL  CCNOITIONS  FOR  TAXI 
C  FuR  j*H  MAIN  SXN  NOSE  GEAR 
C  ZMI=  MAIN  GEAR  TIRE  OEFLECT ION  35  PERCENT 
C  ZNI  =  NOSE  CEAfi  TIRE  OEFLECT 1CN  35  PERCENT  OEFLECTICN 
C  ZCGI=  CG  DEFLECTION 
C  THETAI*  PITCH  ANGLE 

C  XMA IN=  MAIN  LANDING  GEAR  STATIC  STROKE 
C  XNOSE*  NOSE  GEAR  STATIC  STROKE 


RH*H/<l.*A/6) 

RNsM-RM 

RNxRM/SXM 

RN=RN/SXN 

ZHI=-RM/TSM 

ZNI=-RN/T3N 

RSM  =  (RM  -MM)/. 66 

RSN=RN-MN 

XNOSE**PAON* VON/RSN-VON/ AAN 
XMAIN=C*PA0M*V0M/RSM-V0M/AAM-1A.3Z5  1 

IFfPSH.GE. 7669. .ANC.RSH.LT. 45557.)  XHAIN  *  - (RSM ♦2696587.) /189440. 
XMA66=XM AIN/ .66 

TMETAI*- (XNOSE+ZNI-( XHA66+ZMI) ) / (S^A ) 

ZCGI=XHA66-A*THETAUZHI 
00  10  1=1, NFH 
QO(I)=G t 
1C  G( I) =0  . 

RETURN 

ENO 


APPENDIX  II 


FORTRAN  SYMBOL  DEFINITIONS 

This  appendix  ..obtains  an  alphabetical  listing  of  the  Fortran  variables 
used  in  the  piogram  TAXI  categorized  by  the  subroutine  in  which  they  are 
defined.  In  cases  where  a  variable  is  defined  in  two  or  more  subroutines,  it 
is  listed  under  the  subroutine  in  which  it  is  used  most  often.  Some  symbols 
used  in  the  C-5A  computer  code  are  not  listed.  These  variables  are  those  which 
have  been  formed  by  adding  a  1  or  2  to  a  variable  name  which  is  contained  in 
the  basic  1AX1  computer  code.  The  1  refers  to  the  rear  set  of  main  gear  and 
the  2  refers  to  the  front  set  of  main  gear  of  the  C-5A  aircraft.  Thus,  a 
variable  such  as  FSM1  in  the  C-5A  code  .y  be  found  by  looking  for  the  vari¬ 
able  FSM  in  the  listing  of  symbols  and  associating  the  definition  of  the 
variable  with  the  rear  set  of  main  gear  of  the  aircraft.  Some  variables 
denoted  (C-5A  simulation  only)  are  those  contained  in  the  C-5A  computer  code 
exclusively. 
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APPENDIX  II 


TAXI 

SYMBOL 

A 

AAM 

AAN 

AHM 

AHN 

AM 

AMA 

AN 

AREA 

B 

BM 

BhA 

BN 

C 

CD 

CGACC 

CGGUT 

CL 


FORTRAN  SYMBOL  DEFINITIONS 

DEFINITION 

Distance  from  CG  to  rear  main  gear 
Pneumatic  area,  main  gear 
Pneumatic  area,  nose  gear 
Hydraulic  area,  main  gear 
Hydraulic  area,  nose  gear 

Coefficient  of  polynomial  fit  to  runway  profile  segment ; 
rear  main  gear 

Coefficient  of  polynomial  fit  to  runway  profile  segment, 
front  main  gear(C-5A  simulation  only 

Coefficient  of  polynomial  fit  to  runw  ..rofile  segment, 
nose  gear 

Aircraft  wing  area 

Distance  from  CG  to  nose  gear 

Coefficient  of  polynomial  .  to  runway  profile  segment, 
rear  main  gear 

Coefficient  of  polynomial  fit  to  runway  profile  segment, 
front  main  gear  (C-5A  simulation  only) 

Coefficient  of  polynomial  fit  to  runway  profile  segment, 
nose  gear 

Distance  from  CG  to  front  main  gear  (C-5A  simulation  only) 
Coefficient  of  drag 

Array  containing  CG  accelerations  for  Calcomp  plot 
Total  CG  acceleration 


Coefficient  of  lift 


SYMBOL 


DEFINITION 


CM  Coefficient  of  polynomial  fit  to  runway  profile  segment, 

rear  main  gear 


CMA 

CN 

DDPLOT 

DISTAN 

DM 

DMA 

DN 

DX 

DZ 

ELEV 

ELEV1 


Coefiicient  of  polynomial  fit  to  runway  profile  segment, 
front  main  gear  (C-5A  simulation  only) 

uuef ficient  of  polynomial  fit  to  runway  profile  segment, 
nose  gear 

Array  containing  aircraft  distance  down  the  runway  for 
Calcomp  plot 

Distance  down  the  runway  used  in  removing  overall  slope 
from  runway  profile 

Coefficient  of  polynomial  fit  to  runway  profile  segment-, 
rear  main  gear 

Coefficient  of  polynomial  fit  to  runway  profile  segment, 
front  main  gear  (C-5A  simulation  only) 

Coefficient  of  polynomial  fit  to  runway  profile  segment, 
nose  gear 

Time  step  for  Integration 

Variable  which  compensates  for  the  overlap  of  two  adjacent 


Array  containing  runway  profile  elevations 
Elevation  of  first  runway  profile  point 


ENDRUN  Length  of  runway 

GM  Arra*  containing  generalized  masses  for  each  flexible 

mode  of  vibration 

HDR  Counter  for  printing  header  on  printed  output 

[  Index  variable 

IFPLOT  Variable  which  contains  decision  to  produce  plotted  output 

or  i. ot 

II  Subscript  variable  for  runway  markers  on  Calcomp  plot 
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SYMBOL 


DEFINITION 


ITT 

IVAL 

IV ‘1 A 

IXLONG 

J 

LD 

LDI 

LL 

LLL 

LPRIN 

LSD 

LSD1 

M 

MCG 

MGM 

MM 

MMI 

MN 

MRM 

NFM 

NN 


Index  variable 

Integer  truncation  of  distance  between  nose  and  rear  main 
gear 

Integer  truncation  of  distance  between  rear  main  gear  and 
front  main  gear  (C-5A  simulation  nly) 

Integer  truncation  of  length  of  time  axis  on  Calcomp  plot 

Subscript  variable  for  runway  profile 

Counting  variable  for  runway  profile  input 

Counting  variable  for  runway  profile  output  listing 

Subscript  variable  for  storage  of  pilot  station  acceleration 
time  history 

Subscript  variable  for  storage  of  aircraft  speed  and  distance 
for  Calcomp  plot 

Runway  distance  for  runway  profile  listing 

Counting  variable  for  runway  profile  input 

Counting  variable  for  runway  profile  listing 

Counting  variable  for  printing  out  output  header  first  time 

Mass  of  entire  aircraft 

Integer  truncation  of  total  simulation  time 

Mass  of  unsprung  portion  of  one  main  landing  gear 

Pitching  moment  of  inertia  about  aircraft  center  of  gravity 

Mass  of  the  unsprung  portion  of  the  nose  landing  gear 

Length  of  runway  divided  by  1000  feet 

Number  of  flexible  modes 

Subscript  variable  for  CG  acceleration  time  history 
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SYMBOL 


DEFINITION 


NPTSS  Number  of  runway  profile  points 

NSCM  Number  of  slope  or  area  changes  on  main  strut  metering  pin 

NSCN  Number  of  slope  or  area  changes  on  nose  strut  metering  pin 

NTRUN  Defines  run  as  taxi  oi  takeoff 

NiO  index  variable  for  normalization  of  runway  profile 

OAM  Area  of  orifice  hole,  main  gear 

OAN  Area  of  orifice  hole,  nose  gear 

OMEGA  Array  of  flexible  mode  frequencies 

PAOM  Preload  pressure  of  main  gear  strut 

PAON  Preload  pressure  of  nose  gear  strut 

PINDM  Array  containing  main  gear  metering  pin  diameters  for  con¬ 

ventional  aircraft  and  net  orifice  areas  for  aircraft  with 
metering  tubes  or  flulted  metering  pins 

PINDN  Array  containing  nose  gear  metering  pin  diameters  for  con¬ 

ventional  aircraft  and  net  orifice  areas  for  aircraft  with 
metering  tubes  or  flutted  metering  pins 

PLANE  Aircraft  being  simulated 

pr-  Counting  variable  for  storage  ot  distance  and  speed  for 

Cal comp  plot 

PROF  Runway  profile  time  history  elevations 

PROFIO  PROF  (NN+2)X1Q 

PSA  Pilot  station  acceleration 

PSACC  Array  containing  pilot  station  deceleration  time  history 

PSARM  Distance  from  pilot  station  to  CG 

Q  Array  of  non-dimensional  time  dependent  coordinates  which 

weight  the  amount  of  motion  due  to  each  flexible  mode  in 
the  total  motion  of  the  aircraft 


SYMBOL 

DEFINITION 

QD 

Time  derivative  of  Q 

QDD 

Time  derivative  of  QD 

QDDCG 

CG  acceleration  due  to  flexible  motion 

QuDlb 

filot  station  acceleration  due  to  flexible  motion 

QDDTAL 

Tail  station  acceleration  due  to  flexible  motion 

REACTM 

Static,  total  force  at  main  gear 

REACTN 

Static,  total  force  ac  nose  gear 

RM 

Incremented  variable  for  determining  position  of  runway 
markers 

RMARK 

Array  containing  runway  markers  positions 

SICC 

Mo^e  shape  deflection  of  CG 

SI MAIN 

Mode  shape  deflection  at  main  landing  gear 

SI NOSE 

Mode  shape  deflection  at  nose  landing  gear 

SIPS 

Mode  ohape  deflection  at  pilot  station 

SITAIL 

Mode  shape  deflection  at  tail  station 

SITE 

Location  of  runway 

SLM 

Distance  from  CL  of  main  gear  axle  to  CG  of  aircraft 
strut  fully  extended 

with 

SLN 

Distance  from  CL  of  nose  gear  axle  to  CG  of  aircraft 
strut  fully  extended 

with 

SLP 

Overall  slope  of  runway  profile 

SPEED 

Initial  speed  of  aircraft 

SSPLQT 

Array  of  velocity  of  aircraft  for  Calcomp  plot 

STORE1 

Temporary  storage  space  for  CG  accelerations 

SYMBOL 


DEFINITIONS 


STORE2 
STORE3 
STORE4 
S'  SOKM 

STROKN 

SXN 

SXN 

T\ILAC 

TAILRM 

TAKOFF 

THRUST 

TIME 

TIMEX 

TIME1 

TSM 

TSN 

TYPRUN 

VAL 

VALA 


Temporary  storage  space  for  CG  accelerations 

Temporary  storage  space  for  pilot  station  accelerations 

Temporary  storage  space  for  pilot  station  accelerations 

Ar lU'j  of  strut  stroke  values  corresponding  to  metering 
pin  values  (PINDM) ,  main  gear 

Array  of  strut  stroke  values  corresponding  to  metering 
pin  values  (PINDN);  nose  gear 

Number  of  main  gear  struts 

Number  of  nose  gear  struts 

Acceleration  at  tail  station 

Distance  from  tail  station  to  CG 

Rotation  velocity  of  aircraft 

Total  thrust  of  aircraft 

Array  of  simulation  times  at  which  CG  acceleration  time 
history  points  are  stored 

Counter  variable  for  printed  output 

Array  of  simulation  times  at  which  pilot  station  time 
history  are  stored 

Tire  spring  constant,  main  gear 

Tire  spring  constant,  nose  gear 

Define:,  simulation  as  takeoff  or  taxi 

Distance  between  nose  and  rear  main  gear 

Distance  between  rear  main  gear  and  front  main  gear 
(C-5A  simulation  only) 


VOM 


Main  gear  strut  fully  extended  volume 


SYMBOL 

VON 

U 

WM 

WN 

X 

XLONG 

XLONG2 

XPROF 

XSTOP 


DEFINITION 

Nose  gear  strut  fully  extended  volume 
Weight  of  aircraft 
Main  gear  unsprung  weight 
Nose  gear  unsprung  weight 
Simulation  time 

Length  of  time  axis  for  Calcomp  plot 
XLONG/ 2 

Location  for  printing  of  "NOSE  GEAR  TRACK"  on  Calcomp  plot 
XLONG+5 


YP 

YPA 


YPN 

ZADOT 

ZNDOT 


Array  containing  runway  segment  elevation  points  and  slope 
from  end  of  previous  segment,  rear  main  gear 

Array  containing  runway  segment  elevation  points  and  slope 
from  end  of  previous  segment,  front  main  gear  (C-5A  simu¬ 
lation  only) 

Array  containing  runway  segment  elevation  points  and  slope 
from  end  of  previous  segment,  nose  gear 

Slope  of  runway  segment  at  end  point,  front  main  gear 
(C-5A  simulation  only) 

Slope  of  runway  segment  at  end  point,  nose  gear 


I 


li 

il 


TAYLOR 

AOM 

AON 

a)M 

CON 

DRAGA 


Net  orifice  area,  main  gear  (OAM-metering  pin  area) 
Net  orifice  area,  nose  gear  (OAN-metering  pin  area) 
Damping  coefficient,  main  gear 
Damping  coefficient,  nose  gear 
Aerodynamic  drag 
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•SYMBOL 


DEFINITION 


DRAG!' 

FSM 

FSN 

FSTN 

FST1 

FST2 

FTM 

FTN 

F2M1 

F2M2 

F2N 

QTDM 

QTDN 

QTM 

QTN 

SLOPEM 

SLOPEN 

SSM 

SSN 

T(2) 

T(4) 

T(6) 


Rolling  Drag 

Total  force  in  all  main  gear  struts 
Total  force  in  all  nose  gear  struts 

Net  force  on  secondary  piston,  nose  gear  (C-5A  Simulation  only) 

Net  force  on  secondary  piston,  rear  main  gear  (C-5A  Simulation 
only) 

Net  force  on  secondary  piston,  front  main  gear  (C-5A  Simulation 
only) 

Force  in  tires,  main  gear 

Force  in  tires,  nose  gear 

Total  force  in  secondary  chamber,  rear  main  gear  (C-5A  simu¬ 
lation  only) 

Total  force  in  secondary  chamber,  front  main  gear  (C-5A  simu¬ 
lation  only) 

Total  force  in  secondary  chamber,  nose  gear  (C-5A  simulation 
only) 

Total  velocity  due  to  flexible  modes  at  main  gear 

Total  velocity  due  to  flexible  modes  at  nose  gear 

Total  deflection  due  to  flexible  modes  at  main  gear 

Total  deflection  due  to  flexible  modes  at  nose  gear 

Slope  of  line  drawn  through  two  metering  pin  points,  main  gear 

Slope  line  drawn  through  two  metering  pin  points,  nose  gear 

Pneumatic  force,  main  gear 

Pneumatic  force,  nose  gear 

CG  vertical  displacement 

Unsprung  mass  vertical  displacement,  front  main  gear 
Unsprung  mass  vertical  displacement,  nose  gear 
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SYMBOL 


DEFINITION 


T<8) 

Rigid  body  pitch  angle  of  aircraft 

T(10) 

Horizontal  distance  of  airoraft 

T(12) 

Unsprung  mass  vertical  displacement,  front  main 
(C-5A  simulation  only) 

gear 

T(14) 

Vertical  displacement  of  secondary  piston,  nose 
(C-5A  simulation  only) 

gear 

T(16) 

Vertical  displacement  of  secondary  piston,  front 
(C-5A  simulation  only) 

main  gear 

T(18) 

Vertical  displacement  of  secondary  piston,  rear  main  gear 
(C-5A  simulation  only) 

TD(2) 

Time  derivative  of  T(2) 

TD(4) 

Time  derivative  of  T(4) 

TD(6) 

Time  derivative  of  T(6) 

TD(8) 

Time  derivative  of  T(8) 

TD(10) 

Time  derivative  of  T(10) 

TD(12) 

Time  derivative  of  T(12)  (C-5A  simulation  only) 

TD(X4) 

Time  derivative  of  T(14)  (C-5A  simulation  only) 

TD(16) 

Time  derivative  of  T(16)  (C-5A  simulation  only) 

TD(18) 

Time  derivative  of  T(18)  (C-5A  simulation  only) 

TDDC2) 

Time  derivative  of  TD(2) 

TDD(4) 

Time  derivative  of  TD(4) 

TDD(6) 

Time  derivative  of  TD(6) 

TDD(8) 

Time  derivative  of  TD(8) 

TDD(IO) 

Time  derivative  of  TD(10) 

TDD(12) 

Time  derivative  of  TD(12)  (C-5A  simulation  only) 

TDD(14) 

Time  derivative  of  TD(14)  (C-5A  simulation  only) 

TDD(16) 

Time  derivative  of  TD(16)  (C-5A  simulation  orly) 
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SYMBOL 


TDD (18) 

VELM 

VELN 

vli  n 

XMA1N 

XMLK 

XNLK 

XHOSE 

YCEPM 

YCEPN 

Z 

ZPM 

ZPMA 

ZPN 

KM 

RMl'I 

RM2'I 

RN 

RSM 

RSN 

THE1AI 

XMAIN2T 


DEFINITION 

Time  derivative  of  TD(18)  (C-5A  simulation  only) 

Total  strut  velocity,  main  gear 
Total  strut  velocity,  nose  gear 
Aerodynamic  lift  force 
Strut  stroke,  main  gear 
Absolute  value  of  XMAIN 
Absolute  value  of  XNOSE 
Strut  stroke,  main  gear 

Y  intercept  of  line  drawn  through  metering  pin  points, 
main  gear 

Y  intercept  of  line  drawn  through  two  metering  pin  points, 
nose  gear 

Distance  of  aircraft  from  beginning  of  a  4  ft  runway  segment 
Runway  elevation,  rear  main  gear 

Runway  elevation,  front  main  gear  (C-5A  simulation  only) 

Runway  elevation,  nose  gear  (C-5A  simulation  only) 

Static  reaction  force  at  main  gear 

Static  reaction  force  at  rear  main  gear  (C-5A  simulation  only) 
Static  reaction  force  at  front  main  gear  (C-5A  simulation  only) 
Static  reaction  force  at  nose  gear 
RM  -  WM 
RN  -  WN 

Rigid  body  Initial  pitch  angle 

Test  variable  for  rigid  body  initial  conditions 
(C-5A  simulation  only) 


SYMBOL 

ZCGI 

ZMI 

zn: 

COEFF 

A 

.1 

C 

D 

Y 

TLOOK 

I 

N 

P 

S 

SLOPE 

YCEPT 


DEFINITION 

Initial  CG  vertical  displacement 
Initial  tire  deflection,  main  gear 
Initial  tire  deflection,  nose  gear 

Coefficient  of  polynomial  fit  to  runway  profile  segment 

Coefficient  of  polynomial  tit  to  runway  profile  segment 

Coefficient  of  polynomial  fit  to  runway  profile  segment 

Coefficient  of  polynomial  fit  to  runway  profile  segment 

Runway  profile  elevation  values 

Index  variable 

Number  of  values  in  metering  pin  -  stroke  table 

Metering  pin  diameter  or  net  orifice  area  for  aircraft  with 
metering  tubes  or  flutted  metering  pins 

Strut  stroke  in  metering  pin  table 

Slope  of  line  drawn  between  two  metering  pin  points 

Y  intercept  of  line  drawn  between  two  metering  pin  points 
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